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ABSTRACT 


On  August  20,  1978,  the  Syntactic  Aperture  Radar  (SAR) 
or.  board  the  satellite  SEASAT,  gave  3vider.ce  of  high  energy 
internal  oceaiic  activity  at  the  shalf break  in  the  northern 
part    of   the   Bay    of    Biscay. 

Quantitative  spatial  measurements  of  internal  wave 
patterns  were  correlated  with  conventional,  quasi-synoptic 
ir.  situ  da*a  and  yielded  phase  speeds  of  O.SoM/S  for  high 
frequency,  nonlinear  internal  waves,  generated  in  groups  at 
the  canyons  indenting  the  shelfbrea<  and  at  -he  time  of  low 
tide  in  Brest.  Their  oeriods  were  found  to  be  of  70MIN  and 
their  amplitudes  estimated  to  be  15H.  Ocean  swell  refrac- 
tion, observed  on  the  SAR  image,  together  with  a  localized 
significant  increase  in  *ave  height  detected  by  the  SEASAT 
altimeter,  was  related  to  the  shear  of  a  northwestward  gees- 
trophic    jet     (oa.    3.5H/S)     at    "he    break. 

Possible  operational  oonsequenoes  of  these  features  are 
summarized  as  an  assessment  of  a  3AR's  capability  to  give 
tactical  as  well  as  scientific  real-time  information  on  the 
internal    oceai. 


RESUME 

Le      20    Acut      1978,      le      SAR      (Synthetic   Aperture      Radar) 
embarque    sur     le    satellite      SEAS  AT,      a      mis   =r.      evidence   line 
intense   activite      interne   de   1*  ocean      au    voisinage      du   talus 
continental    dans   la    partia    nord   du  Golfe    de   Gascogne. 

Des  mesures  quantitativ as  dss  aanif estat ions  de  surface 
des  ondes  internes  ont  ate  corralaes  avec  des  donnees 
conventionnelles      "in-situ",        et   quasi      simuitanees.  Des 

vitesses  de  phase  de  0,55  metres  par  seconds  ont  ete  mises 
en  evidence  pour  ies  trains  &*  ondss  internes  non-iineairos, 
generes  a  l'haure  de  la  basse  mer  da  Brest,  dans  les  canyons 
bordant  le  plateau  continental;  la  periode  da  ces  ondes  est 
de  70  minutes  et  leur  amplitude  a  ata  estimee  a  15  metres. 
La  refraction  de  la  houle  observaa  sur  1* image  SAR,  en 
conjonction  avec  1* observation  par  l'aitimetre  de  SEASAT 
d*une  amplification  localisee  mais  significative  de  la 
hauteur  de  vague,  a  ete  attribuee  au  courant  de  cisaillemant 
dun  "jet"  geostr ophique  da  surface  (da  l'ordre  la  0,5  metres 
par    seconde)    portant    vers   le  nord-ouest    le    long    du    talus. 

Les  consequences  cperat ionelies  possibles  de  ces  ohenc- 
menes  sont  reca pitulees,  pour  demontrer  les  capacites  d'un 
radar  SAR  a  dDnner  des  informations  tactiques  aussi  bien  que 
scientifiqu as,    en   temps   real,    sur   1' ocean    interne. 


TABLE    OF    CONTESTS 


I.  INTRODUCTION 14 

A.  FOREWARD 14 

B.  DAM     ACQUISITION 14 

C.  OBJECTIVES 16 

D.  UNITS    NOTATIONS    CONVENTIONS       17 

II.  SAR    IMAGING    OF    INTERNAL    WAVES 21 

A.       SYNTHETIC    ANTENNAS 21 

3.       OCEAN    IMAGING    MECHANISMS    2  2 

1.  Doppler   Effects    from    ths    Moving  Ocean 
Surface 23 

2.  Scattering   of   Radar    Signal   from  Ocear. 

Waves 24 

3.  Bragg    Model    versus   Tangent    Plane    Model    .    .  24 

C.  COJPLING    WITH    INTERNAL    OCEAN    DYNAMICS       ....  25 

1.  The    Surface    Straining    and   the    Bunching 
Effects 25 

2.  Application    to    Internal    Waves      26 

D.  SENSITIVITY   OF    \     SAR    TO    INTERNAL    MOTIONS     ...  28 

1.  Environmental  Limitations      28 

2.  Propitious   Internal    tfave    Properties      ...  29 

3.  Favorable    Cases    :     an    Assessment 30 

E.  COMBINED    UTILIZATION    0?    CONVENTIONAL 

SENSORS        31 

1.  Bathythermographs 31 

2.  Current   and    Temperature    Observations    .    .     .  32 

III.  GENERAL    GEOPHYSICAL     FEATURES    OF    THE    BAY    0? 

3ISCAY 35 


A.       ATMOSPHERIC   CIRCULATION      35 

5.        BATHYMETRY 37 

C.  OCEAN    CLIMATOLOGY       38 

D.  WATER    MASSES,     FRONTAL    STRUCTURES,     AND 

EDDIES 39 

S.       TIDES 41 

?.       INTERNAL    WAVE    MOTIONS 42 

IV.  GEOPHYSICAL    SITUATION    ON    AUG05T    20     1978 50 

A.       METEOROLOGY   AND     WIND    STR25S 50 

3.       SWELL 51 

C.  MEAN     AND    TIDAL    CURRENTS 52 

D.  SPATIAL    THERMAL     STRUCTURE 53 

1.  Analysis   cf    the    Hean 54 

2.  Analysis    of    the    Variability 56 

3.  Selected   Representative   Transects      ....  56 
S.       TEMPORAL    THERMAL    ANALYSIS       57 

?.       SUMMARY       59 

/ 

V.  SAR  IMAGE  ANALYSIS 76 

A.   ORBIT  PARAMETERS  AND  GEOGRAPHIC  REFERENCES  .  .  76 

3.   HEURISTIC  ANALYSIS  OF  THE  SAR  IMAGE 77 

1.  Inner  Shelf 77 

2.  Outer  Shelf 78 

3.  Continental  Slope  or  Generation  Area  ...  79 

4.  Continental  Rise  or  Propagation  Area  ...  81 

C.  ANALYSIS  OF  THE  SWELL 82 

D.  SUMMARY 32 

VI.  CORRELATIVE  INTERPRETATION  OF  THE  SAR  IMAGE   ...  39 

A.  TWO  LAYER  MODEL  ON  THE  SHELF .89 

B.  NORMAL  MODES  OF  THE  BAROCLINIC  INTERNAL 

TIDE 91 

C.  HIGH     FREQUENCY    INTERNAL    *AV2S 9U 

1.      Estimation    of  -he    Mixed   Layer   Depth               .  9a 


2.      Estimation   of  the    Internal    Waves 

Amplitude 9  5 

D.       INTERNAL     WAVE    INTERACTION    WITH    SHSLFBREAK 

STUCTURE 97 

S.        SWELL    REFRACTION    AND    SHSLFBREAK    CURRENT       ...  97 

VII.  TENTATIVE    EXTRAPOLATION    TO    3ENERAL    OCEANIC 
FEATURES 104 

A.  INTERNAL     WAVE    FENERATION 10^ 

B.  TIDAL    FIXING    ON    THE    CONTINENTAL    SLOPE        ...  105 

C.  SHELFBREAK   FRONTS    AND    JETS 105 

VIII.  APPLICATION    TO   UNDERWATER    ACOUSTICS       108 

A.       IDEALIZED    WAVE    FIELDS 110 

3.       TRANSMISSION    LOSS    DERIVED     FROM    THE    PS 

MODEL 111 

C.  RAY     TRACING    WITH     MEDUSA 111 

D.  GENERAL    RESULTS ,  113 

IX.  3AR    AS    AN    INTERNAL    OCEAN    SAMPLING    SYSTEM    ....  118 
A.       GENERAL    REQUIREMENTS    FOR    A    3AR    OCEAN 
MONITORING    SYSTEM       118 

3.       IN    SITU    COMPLEMENTARY    SIMULTANEOUS     PROBING    .  120 

C.       OPERATIONAL   PROCESSING    AND    EXPLOITATION       .     .  120 

X.  SUMMARY     AND    RECOMMENDATIONS 123 

BIBLIOGRAPHY      126 

INITIAL    DISTRIBUTION    LIST    132 


LIST    OF    TABLES 

I.  SEASAT    SAR    CHARACTERISTICS    .22 

II.  CLIMATOLOGY    FOR    AUGUST 36 

III.  TIDAL   ZONSTITUENTS .  42 

IV.  METEOROLOGICAL    OBSERVATIONS    3a.47N,005W       50 

V.  COEFFICIENT     AND    TIMES    0?    HI3H    TIDES 5U 

VI.  SURFACE    TIDAL    CURRENTS    AT    47024*  N    006O05'W     ....  55 

VII.  SEASAT    ORBIT     785    PARAMETERS       76 

VIII.  INTERNAL     WAVELENGTHS     93 


LIST    OF    FIGURES 

1.1  GEOGRAPHIC    GRID 19 

1.2  FHOTO    REPRODUCTION    OF    THE    3AR    IMAGE 20 

2.1  SEASAT    SAR    GEOMETRY 33 

2.2  SYNTHETISATION    OF    AN    ANIENS 33 

2.3  DOPPLER    EFFECTS    OF    CURRENT    SHEAR 3U 

2.4  INTERACTION    OF    INTERNAL    AND    SURFACE    WAVES       .     .     .  34 

3.1  TOPOGRAPHY    OF   THE    BAY    OF    BISCAY       ........  44 

3.2  BATHYMETRY    OF   THE    BAY    OF    BISCAY     (SHELF)         ....  45 

3.3  BATHYMETRY    OF   THE    BAY    OF    BISCAY     (OPEN    SEA)      .     .     .  46 

3.4  GENERAL    OCEAN   CIRCULATION 47 

3.5  T-S    DIAGRAM 48 

3.6  COTIDAL    MAP 49 

4.1  SURFACE    ATMOSPHERIC    PRSSSJRE    CHART    60 

4.2  SCATTEROMETER   WINDS 61 

4.3  DAILY    SYNOPTIC    WEATHER    CHARTS       ,     ...  62 

4.4  PROPAGATION    OF    THE    SWELL 63 

4.5  SIGNIFICANT    WAVE    HEIGHTS 64 

4.6  PROGRESSIVE    CURRENT    VECTOR    DIAGRAMS       65 

4.7  DISTRIBUTION    OF    BT    SAMPLINGS     IN    TIME    AND 

SPACE 6  6 

4.3  ALONG    TRACK    AVERAGE    OF    ISOTHERMS     67 

4.9  ALONG    TRACK    STANDARD    DEVIATION    OF    ISOTHERMS       .     .  53 

4.10  OVERALL    TRANSECT    OF    THE    BAY    OF    3ISCAY 69 

4.11  SELECTED    REPRESENTATIVE    TRANSECTS    LOCATIONS       .     .  70 

4.12  SELECTED    REPRESENTATIVE    TRANSECTS       71 

4.13  TYPICAL    THERMAL    STRUCTURE    VARIABILITY       72 

4.14  DYNAMIC    HEIGHTS    AND    GEOSTR0PHY     73 

4.15  DAIL?     AVERAGE    OF    ISOTHERMS    74 


10 


u. 

,  16 

5. 

1 

5. 

2 

5. 

3 

5. 

4 

5. 

5 

6. 

1 

6. 

2 

6. 

3 

6. 

U 

6. 

5 

7. 

1 

8. 

1 

8. 

2 

3. 

3 

8. 

'4 

9. 

1 

I-T    DIAGRAMS 75 

SAR    IMAGE    OF   THE    CONTINENTAL    SHELF 34 

LINE    DRAWING   OF    SURFACE    PATTERNS    85 

SAR    IMAGE    OF   THE    SHELF    BREAK    AND    SLOPE 36 

LINE    DRAWING   OF    SHELFBREAX    AND    SLOPE    AREAS    ...     37 

BATHYMETRY    OF   THE    SHELFBREAS 88 

VERTICAL    MODES    ON    THE    SHELF 99 

VERTICAL    MODES    IN    THE    OPEN    SEA 100 

FIRST     MODES     ADJUSTMENT    TO    D ? EM     SEA    DATA        ...        101 

ESTIMATION    OF   SURFACE    CURRENT    SHEAR       102 

CURRENT    SWELL    REFRACTION    HYPOTHESIS       103 

SCHEMATIC    HYPOTHESIS    OF    SHELF8REAK 

CIRCULATION 107 

SOUND    SPEED    BI-DI MENSIONAL    PROFILES       11U 

HORIZONTAL    ACOUSTIC    PAY    TRACES 115 

TNSONIFICATION  IN  AN  HORIZONTAL  PLANE  ....  116 
TRANSMISSION  LOSSES  IN  AN  HORIZONTAL  PLANE  .  .  117 
OPERATIONAL    USE    OF     SAR    DERIVED    INFORMATION    .     .        122 


1  1 


ACKNOWLEDGMENTS 

I  wish  ro  thank  Professor  !5oo9r=  for  his  having  encour- 
aged ma  to  choose  this  challenging  thesis  topic  and  for  his 
patient      and  dynamic      support      duriig      the    past      year.  He 

frequently  stimulated  new  directions  of  research  and 
suggested  profitable  learning  opportunities;  he  also  intro- 
duced ma  to  many  enriching  personalities  in  the  world  of 
military   and  civilian  oceanography, 

Professor  Hueller  first  inspired  ma  to  work  on  synthetic 
aperture  radar  imagery  and  helped  aa  through  the  difficul- 
ties in  its  interpretation;  his  general  support,  critical 
reviews  and  friendly  encoiragement e  ire  aach  good  reasons 
for    my   gratitude. 

At  the  Fleet  numerical  Oceanography  Center  particular 
thanks  are  owed  to  LCDR  Laath,  LT  Northridge  and  LT  Hergert 
for    both    da^.a  acquisition   and    acoustic   modeling. 

At  the  Jet  Propulsion  Laboratory,  Dr  Le=>  Fu  and  Mr  Ben 
Holt  introduced  ma  to  SAR  data  analysis  techniques;  Mr  David 
Lame  and  Dr  Steven  Pazan  nade  it  possible  for  me  to  use  the 
powerful  Pilot  Ocaan  Data  System  (PDDSI  .  The  visits  I  made 
to    Pasadena    where    both   constructive    and    pleasant. 

I  appreciated  the  opportunities  I  had  to  meet  with  Mr 
Robert  Beal  of  the  Johns  Bopkins  Jiiversity  and  his  helDful 
comments    were   particularly    stimulating. 

Professor  Cantin  pursued  a  tradition  of  generous  supDort 
between  the  Mechanical  Engineering  and  the  Oceanography 
Departments  of  the  School  by  facilitating  access  to  external 
computer    networks. 

I  am  indebted  to  the  French  ^Javy  and  to  ail  the  Naval 
Officers  who  nade  this  challenging  experience  of  workina  at 
the  Naval  Postgraduate  School  possible,  particularly 
Capitaina  de  Corvette  Monsaingeon,  who  was  instrumental  in 
the  approval  of  this  thesis  topic  and  whose  friendly  contact 
was    enriched   by    the    fact    that    he    preceded    me    in    Monterey. 


12 


The  "French        Navy' 5  Service        Hydroarraphigua  -,1- 

Oceancgraphique  contributed  3T  data  and  cliiatDlogic  refer- 
ences in  addition  to  """hose  forwarded  bv  ■the  French  National 
Weather  Service  and  the  U.S.  Navy's  Fleet  Naval  Oceanography 
Center. 

I  would  like  to  express  my  spatial  gratitude  to  Dr 
Raillard  of  the  Centre  National  i' Sxploitat ion  des  Oceans 
for  tha  personal  and  professional  doc a mentation  and  data 
which  he  sent  from  the  "Bangue  des  Donnees  Oceinographiaues" 
which    he    directs. 

This  study  includes  references  from  many  racent  scien- 
tific papers  related  to  the  Bay  of  3iscay  and  to  remote 
sensing  of  internal  waves  in  other  ocean  areas.  Hoping  that 
the  original  intentions  of  the  authors  remain  intact  when 
cited  and  discussed,  I  thank  those  who  sen-1-  reprints  and 
comments  on    these   topics. 

I  dedicate  ohis  work  to  my  Lata  grandfather  Michel 
Homolle. 


13 


I.     INTRODUCTION 

A.  FOREWARD 

Remote  sensing  gives  satellites  an  essential  role  in  the 
monitoring   of  ocean    variability.  Synth-tic    Aperture    Radar 

(SAR)  in  particular,  has  the  ability  ^o  display  images  of 
surface  manifestations  of  internal  processes  which  could 
typically  only  be  inferred  from  laagthy  and  oDstly  in-situ 
sampling    from      ships    at      =  ea   or      instrumented   buoys.  Even 

though  the  dynamical  nonlinear  linkage  involved  in  -he  SAR 
imaging  mechanism  is  presently  incsnplately  understood,  the 
vivid  and  spectacular  images  ;  e.g.,  as  those  acquired  by 
SEASAT,  give  a  tremendous  amount  oi  information  about  well- 
known  or  totally  new  oceanic  phenomena.  To  study  such  an 
image  of  the  3ay  of  Biscay  (Fig.  1.  1  and  1.2)  was  a  chal- 
lenging and  excising  topic  of  thesis  research  for  the  author 
who  has  experienced,  over  several  years,  the  sea-going 
approach  to  military  oceanography  i.i  this  area.  The  excep- 
tional opportunity  for  doing  such  research  was  offered  at 
the  Naval  Postgraduate  School  whers  iirect  support  was  found 
and  the  introduction  to  other  facilities  was  given,  <=spe- 
cially  to  the  Fleet  Numerical  Oceanography  Center  in 
Monterey   and   to    the    Jet    Propulsion   Laboratory    in    Pasadena. 

B.  DATA    ACQUISITION 

The  optically  processed  SAR  image  was  the  starting  point 
of  the  total  data  set  ooilected.  It  was  forwarded  by  the 
Jet  Propulsion  Laboratory  (JPL)  of  the  California  Institute 
of  Technology  in  Pasadena,  and  a  first  heuristic  study  of 
the  imag  =  gave  the  necessary  evidence  of  its  research  possi- 
bilities.       The      swath  is    1  QOKU      wii?    and    -roes      roughly    from 
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Brest  in  Britanny  to  La  Coruna  in  Salicia;  internal  wave 
activity  is  clearly  visible,  particularly  ov=r  the  conti- 
nental slope    (Fig.    1.2). 

Conventional,  quasi-ey noptic  data  were  ^-hen  sought  and 
requested,  either  directly  througa  tha  French  Navy  authori- 
ties, or  by  th21  Naval  Postgraduate  School  from  U.S.  and 
foreign,    military   and  civilian   institutions. 

A      fairly   dense      set    of      XBT's    (bathythermographs)  was 

collected  fr:>m  the  Fleet  Numerical  Oceanography  Center 
(FNOC)  and        the      Franch         "ServLca        Hydr  cgraphique        et 

Oceanographique    da    la   Marine"     (SHOH)  .  The    format    of    these 

BT's  gave  -.he  depth-temperature  pairs  at  turning  points; 
only  five  profiles  extend  below  533M.  Sixty  nine  X3I 
profiles  were  finally  considered,  corresponding  to  a  tima 
period  of  seven  days  around  the  SAR  date  and  located  in  the 
neighbourhood  of  the  image.  The  densest  spacing  covers  the 
continental  slope,  bur  not  where  the  most  intense  activity 
is   noted    from   the   SAR  image. 

Local  atmospheric  da*:a  were  obtained  froi\  the  Fr-nch 
National  Weathar  Service.  Global  *aather  charts  over  tha 
Atlantic  Ocean  and  swell  hindcas*-.  for  the  Bay  of  Biscay  were 
obtained    from   the   FNOC  archive   center. 

Contemporaneous  moored  current  mater  time  series  were 
considered  as  important  auxiliary  information  regardina  the 
internal  baroclinic  modes  of  vertioal  movements,  and  the 
mean  flow.  The  closest  mooring  </as  a  subsurface  array  at 
152U,  3032,  3976J1  depth  anchored  at  46°57«2N,  009057MW 
during  -.he  Northeastern  Atlantic  Daap  Sea  (NFADS)  mission. 
The  digital  data  (current  speed,  iirection,  temperature, 
recorded  depth)  were  forwarded  oy  the  "Centre  National 
d'Exploitation  des  Oceans"  in  Brest,  where  the  "  Bureau 
National   des   Donnees   Oceaniques"    is    locatad. 
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Other  satellite  sensors  carried  by  SEASAT  acquired 
substantial  information  or.  the  sam  =  area  at  rapid  intervals. 
It  was  another  privilege  to  get  direct  computer  access  tc 
the  Pilot  Ocean  Data  System  (PODS)  Df  JPL,  thus  enabling  a 
thorough  survey  of  related  data  given  by  the  3EASAT  altim- 
eter   and   scatt erometer. 

Climatological  data  were  blendel  from  several  sources  in 
the  open  literature,  recent  papers  and  official  unclassified 
documents   of   the    U.S.    and    French    Navies. 

C.       OBJECTIVES 

The  main  aim  in  this  study  is  to  establish  a  definite 
correlation  between  the  SAR  surface  patterns  and  internal 
quasi-synoptic  features  established  by  other  means:  conven- 
tional hydrographic  observations,  current  observations  and 
o^her   satellite    sensors. 

Avoiding  the  complex  problem  of  physical  linkage  in  the 
SAR  imaging  process,  a  pragmatic,  heuristic  approach  to  SAP 
analysis  was  nade,  thus  tentatively  establishing  the  power 
of  SAR  to  monitor  the  oceans.  The  study  has  been  restricted 
to  the  upper  layer  of  the  Bay  of  Biscay  on  each  side  of  the 
continental  slope  and  sheifbreak,  where  both  tha  SAR  image 
and  other  data  bases  were  found  to  provids  rich  information. 
Being  interested  only  in  internal  oceanic  features,  the 
swell  was  only  considered  in  regards  to  interaction  with 
quasi- st eady  currents. 

In  the  coarse  of  the  research,  nany  different  topics  had 
to  be  considered,  thus  greatly  enhancing  the  personal 
benefit  that  could  be  gathered  strictly  from  the  SAR  study. 
Conversely,  this  fact,  added  to  the  limits  in  time  and 
scientific  expertise  of  the  author,  prevents  -his  study  from 
baing      in    any      way    exhaustive.  Firthermore,      an      inhftrent 

feature      of      this      kind      of         quasi-synoptic      study      is      its 
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inability  to  support  objective  "sea-truth.11 .  All  the  results 
presented  her?  must  be  thought  of  as  describing  a  particular 
event  at  a  particular  tine  and  in  a  particular  environment 
that  can  only  be  considered  as  part  of  the  ever-growing 
information   available  to   cceanograpar rs. 

Overall,  the  objective  cf  a  thesis  is  to  increase  the 
knowledge  of  the  author  himself  and,  i.i  this  case,  in  topics 
related  both  to  the  use  cf  conventional  methods  and  of  prom- 
ising satellite-related  techniques.  Consequently,  this 
research  was  a  learning  experience  and  an  enriching  intro- 
duction to  the  new  era  cf  operational  military  satellite 
oceanography. 

D.        UNITS    NOTATIONS    CONVENTIONS 

Unless  explicitly  stated  otherwise,  the  following  units 
will    be   used,    with    the   S.I.    system    as    a    general    basis. 

Units   of   t^me 
GMT    Greenwich  time     (Z)    in   Hours,    Miiates,    Seconds. 
(The    local   summer  time  in   ?rance    is:    GMT-020000; 
the    solar  time    at    the  SAR    image   location    is:    GMT-002000) 
Day-Hour   groups    are    expressed    as    DDHHMM,      implicitly    days    of 
August    1978. 

Units   of_  lenath 
Depths    in   meters    (M) 
Distances    in   Kilometers     (KJ1) 

Units   of  t emoerature. 
Degrees    Celsius     (°C)  . 

Geographic    coordinates 
Geographic   units    are      expressed    in    iscimal    dears-s      of    lati- 
tude   and   longitude     (iat,lon).      Longitudes    are   counted   from   0 
to    360    in    most    confutations. 
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Auxiliary  geographic   gr.id 
An    auxiliary   Cartesian  grii    was      adjusted   over    a    map    projec- 
tion   of      the  area   to      meet    the      needs    of   the      particular    5AR 
image      orientation    and     general      bathymetric   gradient       (Fig. 
1.1)       This    grid      is    in  units   of   kilon eters .  rhe    x-axis    is 

oriented  along  -he  shelfbreak  with  positive  values  towards 
the    southeast    (directed   seaward   in   the   Bay). 

The  y-axis  is  oriented  along  the  center  of  the  SAR 
swath  and  towards  the  northeast  (directed  onto  the  French 
continental   shelf)  . 

Hence  x  values  give  the  relative  positions  on  each  side  of 
the  SAR  image  :(|x|<  60  for  locations  on  the  ground  track:), 
y  values  give  relative  positions  en  each  side  of  the  shelf- 
break. 

The   transformation  was    established   as    follows    : 
Two    geographic    points  clearly   visibls    :n    the    SAR    imaae    and  a 
representative    offshore   bathymetric   sounding      were    chosen   to 
get    the    linear    transform    parameters. 

la-      47.30  x=000.00 

Ion   355.67  y=132.42  Penmarc'h    point. 

la"      U3.59  X=000.00 

Ion  352. 7U        Y=-388.92     Ids  Cairos  point. 

la"   46.00         X=116.7 

Ion   356.00  Y=000.00  Offshore    reference. 

The    corresponding   relationships  ar=    : 

la"    =-0.006261    X    +    0.008075    Y  +    46.73 

Ion   =    0.009205    X    *     0.005620    Y  +    354.92 

X    =-51.315    lat    +    73.73    Ion       -  23772. 

y   =   3  4.05  0    lat    +    5  7.17    Ion       -  24217. 
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Conferral  biconic  projection. 
Circles  1  o  cat  9  XBT"  so  a  ndinars. 
Isobaths    (in   meters)    are    smoothed   to    1QKX, 
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II  •     SAR    IMAGING   OF    INTERNAL    WAVES 

Starting  with  a  survey  of  the  recent  li  tsrature  or.  the 
SAR  principles  and  its  capabilities  -to  detect  oeanic 
features,  the  various  assessments  made  of  the  dominant  mech- 
anisms of  interaction  between  an  internal  wave  field  and  the 
short  wind-driven  capillar y/ultragra vity  waves  that  a  SA? 
can    detect    are    reviewed. 

A.        SYHTHETIC    ANTENNAS 

Synthesizing  a  radar  antenna  consists  of  using  a  moving 
antenna  of  limited  size  and  combining  all  information  gath- 
ered on  a  given  point  of  the  earth  during  a  given  time 
interval  (Fig.  2.2)  .  The  distance  of  travel  on  the  orbit 
gives  the  length  of  a  virtual  antenna  of  much  finer  resolu- 
tion. The  resulting  resolution  is  independent  of  range  and 
speed,  and  it  is  theoretically  half  oh?  real  antenna's  aper- 
ture. In  tha  range  (cross -velocity)  direction,  the  length 
of  coherent  pulses  determine  (in  tas  usual  way)  the  corre- 
sponding   accuracy. 

In  one  case  of  SEASAT,  the  data  were  transmitted  in  real 
time  by  the  satellite  to  one  of  the  =Deciai  processing  earth 
stations.  The  final  record  is  either  a  computer-compatible 
tape  format,  or  a  film  conversion.  The  photographic  output 
is  displayed  in  four  parallel  subswaths  to  account  for  a 
progressive  decrease  in  ohe  radar  response  intensity  with 
range  (Fig.  1.2).  Some  particular  images  over  limited  areas 
have  been  digitally  processed  with  a  valuable  enhancement  of 
re so lut ion . 
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The  extraction  of  a  t wo-dimensional  surface  imag<=  spec- 
trum is  done  either  by  in  optical  Fist  Fourier  Transform 
[Gonzales  1979  and  Shuchman  1978],  or  by  numerical  methods 
[algorithm  of  Hasselman  1990].  Considerations  of  the  rela- 
tive qualities  of  these  processing  techniques  can  be  found 
in    Mitchel    [1982]   and  Tomiyasu   [1978]. 

Characteristics  of  ths  SEASAT  3AR  are  given  here  (Table 
I   and   Fig.       2.  1 )  . 

TABLS    I 
SEASAT    SAR    CHARACTERISTICS 
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B.       OCEAN    IHA3ING    MECHANISMS 


The  application  of  SAR  systems  to  oceanography  requires 
extension  of  the  basic  SAR  earth  imaging  theory  to  account 
for    scattering      from   a     tine- vary in j    surface     that    satisfies 

nonlinear  hydrodynamic  constraints.  On  the  other  hand,  the 
general  sca-tering  models  for  low  resolution  radars  lo  not 
apply      to      the      spac=borne         high-resolution        SAR's.  Th= 

complexity  of  the  problem  results  from  the  widely  nonlinear 
electromagnetic  and  hydrodynamic  interactions.  But  -.his 
nonlinear itv   is      at    the      same    time    the      main    reason      for    the 


remarkable   capabilities      of    SAR      foe    detection      of    mesoscale 
oceanic   phenomena. 

1  .       DoDDlrr    Eg f 9 ens    from   the    Roving    2£±&Z   Surface 

The  fact  that  the  ocean  surface  is  not  rigid  has  two 
effects:  (1|  defocusing  of  the  SAR  image  and  (2)  apparent 
displacement  of  the  radar  echo,  as  cap.  be  seen  by  comparing 
echos  of  moving  ships  and  their  stationnary  wakes  (Cf  para 
V-B-1)  . 

Empirical  studies  have  shown  that  to  :=cov^r  the 
optimal  focus  of  SAR  imagery  over  Dcaan  surfaces,  the  best 
speed  correction  to  use  is  the  long  *ave  phase  velocity,  and 
not  as  would  sesm  intuitively,  the  long  waves  orbital 
velocity  on  which  the  capillary  shocz  waves  ride.  In  the 
case    of    internal   waves,    this   phase    velocity    is    given    by 


<V-9'<^ 


J 


where    »h"   is   the   depth   of   the     thermocline   and   "?4",      "ft."    are 

the  densities  each  side  cf  the  then  ocline.  This  displace- 
ment effect  will  periodica lly  modulate  the  image  intensity 
backscattered  by  long  waves,  causing  in  some  oases  periodic 
dark  and  light  bands  associated  with  the  pattern  of  the  real 
long      wave   oattern.  In   the      rase    of      a    "closing"      current 

(cross- track  component  toward  the  satellite)  "he  surface 
image  will  be  artificially  translated  in  the  positive 
along-track  direction;  hence  if  tiis  "closing"  current  is 
increasing  alDng-track  (fig  2.3),  tie  value  of  this  transla- 
tion will  increase  accordingly  and  there  will  be  stretching 
of    the    image   and    a    lower    3r agg-scat t er    than    expected. 
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2.      Scatter  ing    of  Radar    Sianal   from    3cean    'i.1112 


3ased  on  a  review  [Valenzueia  1978]  of  elsc-rcmag- 
netic  wave  scattering  from  ocea.i  waves,  the  standard 
approach  considers  the  ocean  surface  as  statistically 
described,  bat  a  deterministic  formulation  night  be  mors 
appropriate  when  the  resolution  area  is  small  enough  to 
correspond  to  a  few  wavelengths  oily.  The  former  descrip- 
tion is  nevertheless  easier  and  most  of  the  time  more  accu- 
rate. The  statistical  approaci  nay  lead  either  to  a 
Bragg-resonance-like  mechanism  or  a  "facet"  mechanism 
(tan  gent -plane    theory)    [Brown    1976]. 

3-      Bragg  Ho del    vers  us    Tangent    Plane    Model 

The      3ragg    mechanism      is    11  mi  tad      to    rough      surfaces 
satisfying   the    Rayleigh   criterion    : 


2lcH(rms^C0s8  4 


y 


whare    "H"   is   the   typical   surface      height    fluctua' 
is      the      incident   angle      of      the      radar      wave    or. 
surface.      The  total    baclcsca  tier   is   -hen    given,    by 


ion    and    "v" 
the      wa  ter 


<T°(9».  .  =  Utricle  |g.    (9)|     M>(tCr,«3)    y 


•r      s Dec-. rum    c 


) 
where   "^(Kti,^<j)"    is    the      two-dimensional    ?ouri=! 
the    sea    surface    with    an    ocean      wavenumber    " 2kcos«M    and    "a-  •  " 

represents  *:he  radar  reflactivity  of  the  air-water  inter- 
face. The  alternative  approach  Is  the  tangent  plane  or 
facet  model  applicable  to  gently  indulating  rough  surfaces 
and  small  (§<<  30°)  incidence  aigies  (for  SHASAT  SAR,  9  = 
230)  . 


We    have    in   this    case 
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where  R(0)  is  the  Fresnel  coefficient  for  normal  inciier.ee 
and  "<s2>"  is  the  total  slope  variance  for  ocean  wav^s  of 
long  period.  In  both  cases  ^he  intensity  variations  of  the 
SAR  signal  will  be  corr-iated  to  the  energy  density  ,'H,"  of 
the  ocean  waves.  (They  will  be  at  the  same  time  modified  by 
the  mean  incidence  angle,  resulting  in  a  decreasing  average 
intensity  along  a  cross  track  line.)  Therefore,  the  nonli- 
near periodic  energy  tranfer  from  ths  capillary  wav=s  to  ^he 
mesoscale  phenomena  (currents,  swell  and  internal  waves) 
will    dramatically  affect    the    SAR    iaagery. 

C.       COUPLING    WITH    INTERNAL    OCEAN    DTSAHICS 

Small  coherent  disturbances  in  t.ie  current  field  can 
cause  important  modulations  in  the  ultragra vity  wavefield 
[Hollo-Christensen  1982]  :  e.g.,  for  wavelengths  of  30  cms, 
a  shear  current  Df  0.003  3/sec  can  induce  a  modulation  at 
the    surface    to    which    a  SAB     will    be    sensitive. 

1  •      The    Sur f ace    Strain! ng   and   the    3unchina   Effects 

The  possible  mechanism  involved  appear  to  be  either 
an  hydro dynamic  effect  (the  surfacs  straining),  or  a  Doppler 
effect    (the    bunching). 

The  straining  effect  applies  to  all  infernal  ocean 
phenomena  (internal  waves,  shear  currents,  eddies). 
Variations  of  "  0"  n  along  the  direr-ion  of  wave  propagation 
are  induced  by  hydrodynamic  interactions  between  short  reso- 
nant waves  (A  -  • 3M)  and  large  ocean  periodic  motions 
(A>100M).  This  mechanism  can  be  thought  of  as  an  alter- 
nating pattern  of  divergence-convergence  of  the  free 
surface,  riding  on  a  long-wave  modulation  (swell  or  internal 
wave)  .  The  resonant  waves  are  saturated  locally  at  the 
region  of  maximum  convergence  of  velocity,  aheai  of  the  wave 
crest    (Fig.    2. 4)  . 
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In  the  short  wave  (Phillipsi  spectrum,  the  modula- 
tion transfer  function  oan  be  written  [&lpersrRoss  and 
Ruff enach   1981  ]: 


V=Asil 


M>o  cb 
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where  "H"  is  a  time  constant  relatad  to  the  resonant  waves, 
"<p"  is  the  angle  between  the  flight  (azimuth)  direction  and 
ths  direction  of  propagation  of  the  internal  wave,  and  "JL" 
is    the   angular    frequency    of    the    wavs    component. 

The  vslocity  bunching  effect  is  the  result  of  arti- 
ficial displacement  of  the  image  of  periodically  moving 
surface  [Aipers  and  Rufenach  1979,  Mpers  et  ai.  1981]. 
Simply  put,  the  Doppier  shift  die  to  surface  velocity  is 
attributed  to  location  (azimuth)  and  the  siar.al  is  mislo- 
cated    in      the   image.  It    seems      to    be      the    best      model   for 

azimuth  traveling  swell  and  eventually  for  internal  waves. 
However,  it  is  a  geometric  effect  that  doss  not  account  for 
^nsrcy  transfer  as  wcll  as  ~he  surface  straining  dynamic 
mechanism. 

2.      Application    to   Internal   Waves 


»tfa" 


or      a    wave      component    or 


The    general      form   o: 
wave    number   "KM ,      angular   frequency    "&" ,      and   amplitude   "A", 
may      therefore    be      written    more      generally      as      Vesecky      and 


Stewart    198  2  1: 


<T3=  <  <r°>[  1+  JL  cos(kT-H+^)] 
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with 


i  =    arctan  | 
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where    "B"      is   the   total      modulational    transfer      function    and 
can    be   considered   as    the      result    of    two    different    mechanisms 
(neglecting   the    tilting   effect   applicable   only   to    swell) : 
R    =      Rfhydrodyn.)     *    R  (vel.  bunching) 
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With  a  theory  relating  the  surface  wave  field  to  the 
internal  wave  amplitude  and  surface  winds,  information  on 
the  internal  «av?  amplitude,  and  simultaneously  on  the  wind 
field,  can  be  obtained  from  the  imagery.  However,  such  a 
theory  is  not   available   yet    [3eal    1931]. 

It  is  generally  believed  that  internal  waves  are 
seen    by    a   radar    as    a    result    of   two    nechanisms    : 

(1)  The  surface  currents  associated  with  the  internal  wave 
field  sweep  together  surface  oils  and  materials  in  the 
surface  water  convergence  zones  resulting  in  severe  damping 
of   short    capillary/gravity    waves    to    form    slicks,    and 

(2)  The  surface  current  field  induce!  by  -.he  internal  waves 
result  in  a  periodic  modula-ion  of  the  height  of  the 
capillary/gravity   waves. 

actually  the  two  possible  mechanisms  implied  give 
contradictory  results:  damping    3y      oil    slicks      will    occur 

over  the  convergence  lines,  whereas  the  current  hydrodynamic 
damping  will  be  localized  on  the  divergence  areas.  It  is 
the  author's  assamDtion  that  whene/er  the  surface  patterns 
visible  en  the  S  AR  image  are  linear,  localized  and  filamen- 
tary (Cf  patterns  describe!  in  7-B-ll  ,  the  oil  slick  expla- 
nation is  appropriate,  whereas  alternation  of  high  and  low 
brightness  can  be  mere  easily  explained  by  hydrodyr. amical 
damping  due  to  internal  waves  (Cf  group  pat-ems  described 
in    V-B-3    and   7-B-U)  . 

According  to  the  surface  straining  mechanism,  an 
internal  wavs  traveling  at  phas=  speed  " c"  produces  a 
surface  straining  and  a  velocity  pattern  '•u(x-ct)  "  moving 
with    the    waves.  Short    surface   gravity    waves     traveling   in 

the  same  general  direction  as  the  internal  waves  ear- 
strongly  interact,  particularly  when  the  group  vslociny  of 
the      surface  wave      is  close      no      ths    ohas=      velocity    of      the 


internal    wave. 
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Gargett    and      Hughes    [1981]  give      another    description 

cf  the  sani=  principle.  The  horizontal  surface  currents 
associa-ed  with  internal  waves  are  Large,  being  of  the  same 
order  as  the  wave  phase  velocities  themselves.  They  are 
oriented  approximately  either  parallel  (in  the  troughs) 
antiparaliel  (on  the  crests),  to  the  direction  of  propaga- 
tion. A  wind  field ' superi mpcsed  on  tie  current  will  excite 
a  capillary/ultr agravity  wave  spectrum  that  is  steeper  where 
the  wind  opposes  the  internally  generated  current  than  where 
it  is  parallel  to  it  (Fig.  2.U).  Eventually,  a  modulational 
interaction  couples  the  surface-wave  spectrum  to  that  of  the 
internal  waves  and  possibly  creates  a  growth  of  the  latter 
[Dysthe    1981]. 

Indirect  detection  of  internal  waves  may  be  possible 
by  analyzing  the  variations  of  ths  swell  energy  spectrum 
directly  extracted  by  FFT  of  the  SAR  data.  Refraction  of  the 
swell  is  then  the  mechanism  by  whioh  tha  internal  waves  are 
manifested  on   the   sea   surface. 

D.       SENSITIVITY    OF    A    SAR    10    INTERNAL    MOTIONS 
1.      Environmental  Limitations 


For  possible  raiar  imaging  of  ocean  internal  wave 
motion,  the  first  essential  condition  is  the  existence  of 
surface  capillary  or  ultra -cravity  waves.  This  excludes 
total  calm  wind  conditions,  but  not  rough  and  even  breaking 
seas.  The  relative  directions  of  the  short  surface  waves 
and  the  internal  wave  field  affects  the  hydrodynamical 
interaction  process.  Furthermore,  the  displacement  of  the 
satellite  may  induce  DopDler  imaging  effects  dependent  on 
the  angle  which  the  radar  waves  make  with  the  surface  waves. 
Hence  ,  for  a  given  internal  wave  field  and  known  satellite 
orbi-  parameters,  the  optimal  wind  velocity  conditions  could 
be    stated.      It    seems,    but    it    would   need    experimental   confir- 


nation,  that  a  wind  blowing  parallel  to  the  propagation  of 
the  internal  waves  is  optimal  for  strong  interaction.  The 
radar  modulation  is  largest  for  internal  waves  traveling  in 
azimuthal  direction,  but  on  the  other  hand,  it  is  assumed 
that  a  satellite  track  parallel  to  the  same  internal  wave 
field    may   indice   enhancing    Doppler-ef f ects . 

That  these  geophysical  factors  are  sometimes  contra- 
dictory does  not  seem,  however,  to  inhibit  the  capability  of 
a  SAR  to  display  surface  manifestations  (observed  in  many 
different  casss)  of  internal  motions.  In  particular  the 
analysis  conducted  in  chapter  V  shows  internal  waves  propa- 
gating in  almost  all  relative  directions,  thus  proving  that 
at  least  some  interaction  and  radar  modulation  exists  for 
every   ccnfiaurat ion. 

2.      Pro pit i ous    Inter nal    £ave    ProDerties 

The     group    velocity      of   wind- driven      short    waves     is 


given    by 


2  .   J-  +  kT 


3  "  Tic         V 

where   "T"   is   the   sir  face    tension.      The   optimal   internal   wave 
phase   speed    should    then    bs    of    the    same    order    such    that 

implying    approximately   for    ult ra-grav in y   waves 

M   ^      >      j 

which,      in    ths    case    of  the    30CM   wavelength    waves    that   SEASAI 
SAR    can    detect    gives 

AG     ^       i 

— — —     f\»       

For    a   two-layer    ocean   with    a   thermocline    an    5  OS    we    obtain 

AL  ^  o.s  \o~\ 


The  surface  horizontal  particle  velocities  mist,  or.  the 
other  hand,  be  large  enough.  But  *s  know  they  are  directly 
related   to   the    amplitude    of    the   intarnal    wave    itself. 

We  therefore   think    that      detection   of   internal    waves 
by    their    surface   manifestation   should    be    favored    by: 
a    deep   mixed  layer    (H), 
a   sharp   thermocline,    and 
a    density    jump   of    the    oris;:    of    1    /   4  0    H 
(H    in   meters) 

3.      Favorable   Ca  ses    :    an    As3f_ssa?nt 

Internal  tidally  generated  waves  on  -he  continental 
shelf  have  frequency  properties  aad  amplitudes  thai  cause 
them  to  fall  into  the  stated  favorable  conditions.  Those 
internal  waves  are  generated  most  oo urn  only  in  the  highly 
stratified  seasonal  thermocline  by  interaction  between  a 
barotropic   internal    tide    and   bathyuatric   irregularities. 

As  a  particular  example,  the  continental  shelf 
offshore  France  in  the  3ay  of  Biscay  has  strong  tidal 
currents  relative  to' a  semidiurnal  tide  propagating  north- 
ward. The  day  (20  August  1978)  Df  passage  of  SEASAT  over 
■•This  area  corresponds  to  a  spring  til*,  and  the  bathythermo- 
graphs show  the  strong  stra tif ication  oommon  to  the  summer 
season.  A  search  of  meteorological  data  has  confirmed 
average  summer  conditions  to  have  pravaiied,  with  probable 
winds   from      the    northeast.  Previous    conventional      studies 

have  shown  that  internal  waves  ar=  frequently  observed  and 
propagate  perpendicular  to  the  Vlif/SE  orientation  of  the 
continental   slope. 
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E.        COMBINED   (JTILIZATION    OF    CONVENTIONAL    SENSORS 

Adequate  knowledge  of  wind  conditions  is  necessary  for 
tha  reliability  of  any  s AR  analysis.  Ship  observations,  is 
well  as  data  from  other  rsmot?  sensors  (s.g.,  a  scatterom- 
eter)  ,  were  used  in  this  study.  (Both  the  scatt erometer  and 
SAB  are  side-looking  sensors.  Henoe  they  offer  generally 
overlapping  coverage.)  The  required  knowledge  of  -he  stra- 
tification and  of  the  spectrum  of  periodic  movements  can  be 
derived  from  bathythermographs,  and  from  moored  instrumenta- 
tion   measuring    current   and    tempera-are   time-series- 

1  •      Bathythermographs 

One  bat hythermcgraphic  or  ba t hycelerimet r ic  measure- 
ment allows  the  computation  of  tha  eigenvalues  solution  to 
the  internal  *ave  generation  problem  * Baines  1982].  It  is 
then  possible  to  find  the  wavelength  of  internal  waves  from 
the  spacing  of  wave-packets  generated  by  tidal  flow  over 
topography. 

If  t  hs  repetition  rate  of  the  forcing  "r"  is  known 
(tidal  internal  waves  for  instance)  one  can  derive  from  the 
spacing    between    packets    "d" ,    the    piiss    speed 

Cf  -  cl  r 

A  bathythermograph    in  situ    makes    possible   the    computation   of 
the    main    density   gradient    on   the   thsrmocline    from 

4  =  3*  f   • 

Comparing      observed      and      theoretical    values      of      the      phase 
speed,    it    is  possible   to   infer    the   thermocline   depth    by 


u  - 


X2r*        P 


Ap 


/ 
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and  e ver.tuaiiy ,  by  combining  this  information  with  sea- 
surface  temperatures,  tha  heat  storage  in  the  ocean 
[ Hollc-christ=ns en    and  Mascarenhas   1979]. 

2-      Current    and    Tempera  ture   Observations 

Time  series  of  current  and/or  temperature  measure- 
ments at  several  depths  make  possible  the  vertical  modal 
analysis      of  the      vertical      and      horizontal   movements      at      a 


poi: 


•s  nt 


Hiah      resolution      o bservation s    of      the      sea      surface 


temperature  by  airborne  or  spaceborie  infrared  radiometers 
might  detect  horizontal  gradients  dae  to  variations  of  the 
surface  thermal  layer  depth  directly  related  to  the 
stretching  effect        of  an  internal        wave  field. 

Quasi-simultaneous  temperature  transects,  in  the  sama  way, 
can  resolve  the  space  characteristics  of  an  internal  wave 
field  that  can  only  be  inferred  from  a  fixed  observation 
station    together   with  internal    wav?    theory. 


3  2 


^Spacecraft  velocity 
'    vector  v  —  8  km/s 


Synthetic  aperture, 
length 

L  ~  16  km 


Nadir 


D,  real  aperture,  —  12  m 


i9  =  real  beamwidth 
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=  one-half  real  aperture 


From    Scacaborr.e   5A3    £:r   Dcsaaagraphy   [1981] 
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Figure   2.3         D0PPL2R    3FF3CTS    3?    CUHSSST    5H3A3. 


Dashed  linss    are    isotachs 

Snail   arrows    indicate    water    aovea^nt. 
"Co"    is   tha    phase    sos^i     (orooaaation 
Black    arrow    points   to    surface    fio. 
3hioe   arrow    points  to    surface    iaipi.ig 


f    the   isothertnsl 


Figure    2.  '-i         I  !TI  2  3  ACTION    OF    INT3HHAL    HID    SIJ3F\C2    "?i72S. 
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Ill-     GENERAL    GEOPHYSICAL    FEATURES    OF    THE    BAY    OF    BISCAY 

According     to      Sverdrup      et   al.         J1942]      "the      terminal 

branches  of  tie  Gulf  Straain  system  ire  not  all  well  known; 
irregular  branches  which  turn  to  the  south  terminate  in 
great  whirls  off  the  European  coast".  These  irreaularit ies 
could  correspond  to  semi-permanent  circulation  patterns  or 
propagating  eddies.  Recent  studies  have  be  =  n  made  of  the 
mean  circulation  and  of  mesoscale  features  in  the  Bay  of 
Biscay,  particularly  in  the  vicinity  of  the  Armorican  conti- 
nental slcpe,  from  which  some  general  conclusions  can  be 
dr  a  w  n . 

Several  general  geophysical  features  relevant  to  the 
subarea  covered  by  the  SAB  image  will  be  considered,  giving 
particular  emphasis  to  the  summer  cliatatological  conditions 
under    which    the    SAR    image    was    taken. 

A.       ATMOSPHERIC    CIRCULATION 

Atmospheric  conditions  over  the  3ay  of  Biscay  are  deter- 
mined by  the  Location  of  the  North  Atlantic  Subtropical  High 
(pressure  center)  centered  in  February  around  25°N  and  in 
August  around  35°N.  The  High  extends  over  20  to  25  degrees 
of  latitude  in  summer  (vs.  only  15  in  winter)  and  is  limited 
to  the  north  by  the  Icelandic  Low.  The  general  atmospheric 
circulation  ie  consequently  a  mean  eastward  flow,  weak  in 
summer,  which  is  modified  by  the  presence  of  the  European 
continent  and  perturbed  by  the  paeeage  of  fronts  originated 
over  the  Nortwest  Atlantic  in  the  vicinity  of  the  Polar 
Front  (between  45°and  50°M  in  summer).  Table  II  aives  the 
s-atistical  averages  for  August  at  three  stations,  resD<=c- 
fively  north,    east,    and    southeast    06    the    Bay    of    3iscay. 
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TABLE-  II 
CLIMATOLOGY  FOR  &(JSu"ST 
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B.       BATHYMETRT 

The  Eay  cf  Biscay  (Golfs  de  Gascogne)  is  a  fairly  large 
bay.  It  is  450KM  (250NM)  wide  at  its  western  opening  or.  the 
West  European  Basin  and  extends  553KM  (300  NM)  to  the  south- 
east (Fig.  3.1).  Its  physiographic  characteristics  hav- 
been    presented    by  Laughton    [1975]. 

The  abyssal  plain  has  an  average  depth  of  4500M, 
increasing  southward  to  490 0M.  It  is  limited  to  the  south 
by  the  narrow  (40KM)  Cantabrian  continental  shelf,  and  ^o 
the  north  by  a  wide  (200KM)  and  shallow  (130M)  shelf  (Fig. 
3.2)  .  This  latter      shelf    extends      far    beyor.i      thc    Bay      of 

Biscay  to  encompass  the  British  Isles.  The  designation  of 
Celtic  Sea  is  generally  giver,  to  the  northern  parr  of  the 
shelf;  the  Armorican  Continental  Shelf  is  to  the  south  of 
48<>N,    the    latitude    of   Brest. 

The  "La  Chapelle  Bank",  at  tie  above  approxima-e  limit 
(48QN,  007QW)  is  100  to  150  M  deep  aid  characterized  by  large 
scale  and  smooth  ridges  (amplitude  oa.  40M)  oriented  NNF-SSW 
and  20K1  apart.  Its  southwestard  Liiit  (47041 »N,  007013»W), 
sometimes  referred  to  as  the  bank  per  se,  exhibits  sand 
waves  [Cartwright  1958],  correlated  with  scattering  layer 
undulations  [Stride  1960]  is  a  consequence  of  internal  tidal 
wave  activity  " Carruthers  1963  and  Cartwright  1963].  These 
sand-waves  have  an  average  height  of  8M  and  a  wavelength  of 
850H;  •'■.hey  are  oriented  toward  110°T  (parallel  to  the  edge 
of  the  shelf)  and  extend  over  at  least  20KM  from  the  slope. 
Southeast  of  La  Chapelle  Bank,  the  continental  seafloor  has 
a  variable  geological  nature  (rocks,  mud,  gravel  and  sand)  , 
homogeneous    in    depth    with   2.    v=ry    gentle    slope    (0.05%). 

The  shelfbreak  of  the  Armoricar.  Shelf  is  oriented  130°T 
and  occurs  between  180  and  200M  (Fig.  3.3)  The  width  of  the 
continental  slope  decreases  from  35K!1  at  La  Chapelle  Bank  to 
20KM    at    470N,    3  05O30»W,       where    the    steepness    is    maximum    (20*, 
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between  the  233  and  1000H  isobaths')  .  This  steep  slcp=  is 
indented  by  many  canyons  and  ridges  normal  to  its  axis.  The 
designations  given  by  Laughter  [1975]  will  be  ased  here.  On 
the  swa-.h  covarsi  by  the  SAR  image,  five  major  canyons  ara 
visible:  the  Penmarc'h  Canyon,  :is  Suilvinec  Canyon,  the 
Odet  Canyon,  the  Blavet  Canyon  and  the  Be  lie* I  Is  Canyon 
(Fig.    3.  3)  . 

The  continental  rise,  south  of  the  slope,  is  200KM  wide. 
It  is  broken  by  the  Gascc?.y  Knoll  at  45O20'N,  005°20,W,  with 
a  minimum  depth  of  3950.1,  at  ths  =astsrn  adgs  of  the  SAR 
image.  The  3a.ntabr  iar.  Knoll  centarsi  at  450Q5'N,  OOBOOO'W 
is  steeper  and  higher  (3003 M  minimum  iapth,  20  by  40KM  hori- 
zontally) ,  bat  it  is  50KM  off  ths  western  sdge  of  the  SAR 
image. 

C.       OCEAN   CLIMATOLOGY 

The  surface  circulation  in  ths  Northeast am  Atlantic  is 
not  as  well  defined  as  that  in  :i=  western  (Gulf  Stream) 
regime  of  the  North  Atlantic  anticyrion ic  circulation.  After 
crossing  the  mid-Atlantic  ridge,  :i9  Gulf  Stream  lesss  ins 
coherency  and  splits  into  several  Dra.ichas,  which  are  weak 
and  greatly  variable  en  all  space  and  time  scales  and  which 
eventually  decay  into  fields  of  poorly  defined  eddies  (Fig. 
3.4) .  The  easternmost  of  these  branches  flows  south  along 
the  European  roast,  probably  creati.ig  a  weak  cyclonic  circu- 
lation in  the  Bay  of  Biscay,  as  it  waakens  along  the  north 
coast  of  Spall  [Millet  1930].  A  transient  westward  surface 
current  in  ths  northern  part  of  ths  Bay  follows  long  periods 
of  eastward  winds;  i*  is  referred  to  as  the  Ren n el  Current 
and  described  in  the  French  "Instructions  tfautiques"  [SHOM 
1975]   as    being    probably   narrow,    but    strong     (0.5    to    0.73/S). 
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The  near-bottom  circulation  is  highly  variable,  too. 
However,  ^-.here  is  a  quasi  permanent  cyclonic  bottom  circula- 
tion on  the  abyssal  plain  centered  at  44O30'N,  004O30'W, 
which  is  driven  by  the  eastward  flow  along  the  Iberian 
Peninsula.  3n  the  Armorican  Shelf  the  mean  bottom  current 
drifts  to  the  northwest  from  March  to  July,  and  then  towards 
the  southeast  from  August  to  November.  The  average  current- 
speed  measured  at  1M  above  the  sea  floor  is  0.04M/S.  [Camus 
1971  ]. 

The  intermediate  waters  in  the  abyssal  plain  of  the  3ay 
are  under  the  influence  of  the  Mediterranean  outflow,  which 
spreads      north      out      cf   the      Straite      of      Gibraltar.  This 

intruding  water  mass  is  o ante  red  it  a  depth  of  abcu-  1100M 
and  its  T-S  characteristics,  origiiaLly  1=11.9,  s=36.50°/oo, 
decrease  by  mixing  to  T=19,  S=35.3D°/oo  in  the  greater  part 
of  the  Bay  of  Biscay.  Assuming  a  level  of  no-motion  at 
2009M,  results  from  two  separate  hydrographic  data  sets 
(July  196  3  and  November  195  7)  exhibit  the  following  charac- 
teristics of  the  Mediterranean  Water  circulation  [ Le  Floch 
1969]:  a    permanent     cyclonic      circulation      occurs   in      the 

southwest  portion  cf  the  Bay  cf  Biscay,  approximately 
centered  at  4403Q'N,  G07O30'W  on  which  is  superimposed  vari- 
able meanders  and  probable  eddies;  the  speed  of  this  circu- 
lation   is   of  the   order  of   9.05M/S. 

D.        WATER    MASSES,    FRONTAL    STRUCTURES,     AND    EDDIES 

Hydrograph ically ,  the  deep  ocean  waters  of  the  Bay  of 
Biscay  have  three  main  components;  the  relatively  warm 
surface  waters  are  from  the  North  Atlantic  ant  icy clonic 
gyre;  the  deep  and  bottoai  waters  are  derived  from  Antarctic 
Bottom  and  Intermediate  tfa-.er;  tie  midcolumn  water  mass 
char cteri sties  vary  with  the  degree  of  mixing  of 
Mediterranean      Waters     with      either    one      of      the      first      two 
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components.  On  the  Arm  or  i  car.  Shelf  a  mass  of  colder,  less 
saline      water   occurs      as    a  bottom    layer      adjacent    _o      the 

French  coast    (Fig.    3.5). 

we  will  CDnsider  as  a  front,  the  locus  of  sharp  hori- 
zontal gradients  in  temperature  or  salinity  at  an  arbitrary 
depth.  Two      main      surface      fronts      are      driven      by      tidal 

currents,  one  offshore  Britanny  due  to  a  horizontal  decrease 
in  mixing,  and  the  other  at  the  continental  slops  which  will 
be    discussed  in    mors    detail.  Two    pelagic    fronts    have   been 

identified      :  (1)the      the rmalhalii s      front    delimiting      the 

extension  of  Mediterranean  Slaters,  which  is  difficult  to 
monitor  but  which  may  havs  important  acoustical  effects,  and 
(2)  the  thermal  front  bounding  the  ooastai  bottom  cold  water 
mass  of  the  Armorican  Shelf,  which  is  enhanced  in  summer  by 
surface  heating,  but  which  has  no  surface  temperature  mani- 
festation  in  that   season. 

Eddies  occur  in  the  Bay  of  Bisoay  in  summer  conditions, 
when  the  thermal  structure  of  ths  upper  water  column  is 
strongly  developed.  Basel  on  infrared  satellite  images  of 
eddies  bordering  ths  Ushaat  frontal  zone,  these  eddies  are 
most  apparent  along  geostrophically  balanced  fronts  where 
cross  frontal  mixing  takes  the  form  of  cyclonic  disturbances 
[Pingree  1978].  They  have  a  length  scale  of  20  to  no  KM, 
with  a  time  scale  of  a  fsw  days  and  translational  velocity 
less  than  2CM/S.  In  the  case  of  the  Armorican  she  If break, 
irregular  eddies  form  on  each  side  of  ths  200M  isobath;  in 
this  process  "the  cooler  water  sinks  below  the  warmer  water 
and    is   stretched   vertically"   [Pingree    1978]. 

Over  the  abyssal  plain,  less  apparent  eddiss  have  been 
detected  by  drif ting-buoy  experiments  and  by  satellite 
imagery  (Diskson  1931).  These  eddies,  both  anticyclonic  and 
cyclonic,  seem  to  partly  decouple  salinity  and  temperature 
features  from  density  surfaces.  From  an  AVHRR  imag=  taken 
on    14    May    1979,         and  in    the   vicinity    of    the      SAP.    imaae,      it 
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appears  that  cyclonic  eddies  were  aligned  en  tha  continental 
rise  east  of  006°W,  and  that  a  pale  of  cyclonic  and  ar.ticy- 
clonic  eddies  was  located,  resoecti vely ,  northwest  and 
southeast  of  45O30'N,  007O30,».  Tie  length  scale  cf  these 
deep-water    eddies   is    50    to    100KM   [Dickson    198  1], 


E.       TIDES 

"The  tides  of  the  Atlantic  Ocean  can.  be  considered  as 
composed  of  two  parts,  a  co-csciliat ing  tide  which  is  main- 
tained by  the  tida  of  the  Antarctic  Ocean,  and  a  free  fide 
which  is  maintained  by  the  direct  effect  of  the  tide 
producing  forces.  As  a  result,  tha  semidiurnal  tide  can  be 
considered  as  composed  of  two  standing  oscillations  having  a 
phase  difference  of  one  guarter  of  three  lunar  hours." 
[Sverdrup  1942]  The  amphidromic  point  is  "estimated  reli- 
ably" [Cartwright  1930]  at  5QON,  039ow,  so  that  the  tida 
will  propagata  from  south  to  north  la  the  Bay  of  3iscay  with 
a  90  degrees  Greenwich  phase  lag  in  the  south  versus  110  in 
the  north  (Fig.  3.6).  The  tidal  amplitudes  increase  from 
1.0M  offshore,  to  1 . 2M  over  the  continental  rise  and  1.4M  at 
SOK^I  off  the  French  coast.  By  way  of  comparison,  the  diurnal 
tidal  ampli-.udes  over  ths  area  are  of  the  order  of  0.07m 
(Table  III) .  The  tidal  energy  dissipation  computed  by 
Cartwright  [1930  1  largly  exceeds  direct  estimates  based  on 
bottom  friction,  leading  lim  to  conclude  that  "a  substantial 
loss  must  occur  along  the  Biscay  shelf  south  of  Ouessanf 
(Ushant  Island)."  And  indeed,  pingrae  "  198  2]  infers  from 
both  observations  and  nonlinear  tidal  models  -hat  "aboa*  60% 
of  the  tidal  energy  flax  onto  ths  shelf  occurs  in  the 
regions  where  the  tidal  streams  reach  about  49CJVS  (ca. 
005°30'W    to      303 033 'W)".  Furthermore,       he   notes      that    "the 

shelfbreak   bathymetric  contours   tend      to    increas*    this    tidal 
strength    wherever    they  ace    concave    (canyons)". 
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Internal  tidal  modes  have  been  studied  from  current 
meter  observations  at  47Q47N  010<M3'»  in  April  1978,  [James 
1982].  Results  support  the  notion  that  internal  tides  are 
dominated  by  the  lowest  vertical  nodes  and  that  the  M2 
component  is  significantly  coherent  with  t ha  equilibrium 
tide.  Barotr^pic  tidal  current  ellipses  war?  found  to  ba 
narrow  and  directed  perpendicular  to  the  shelf,  thus 
"forcing  the  isopycnal  surfaces  3/sr  topographic  features 
(the  accepted  explanation  for  the  generation  of  internal 
tiles)"  [James  1932].  They  show  a  regular  turning  (counter- 
clockwise) of  their  orientation  axis  from  570  to  3670f!  depth 
and  could  thus  produce  baroclinic  activity.  The  tidal 
ellipses  are  themselves  of  clockwise  rota-ion  at  depths  less 
than    1270M. 


F.       INTERNAL   »A7E    MOTIONS 

Internal  waves  are  a  oommon  phenomenon  in  the  stratified 
ocean  and  exist  in  the  frequency  band  from  the  local  iner- 
tial  frequency  (16. 4H  pariod  at  as^Mi  to  the  maximum  Brunt- 
Vaisala  frequency.  At  tidal  frequencies  (e.g.,  such  as, 
semidiurnal),         these      oeriodic      motions      should      be      called 
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internal  (baroclinic)  tiles  as  they  result  from  forcing  by 
tidal  forces.  At  frequencies  greater  than  tha  Brunt-Vaisala 
frequency,  perturbations  cannot  propagate  as  waves  and  are 
turbulent. 

In  the  canter  of  the  Bay  of  Biscay,  temperature  time 
series  [  Hecht  and  Hughes,  1971]  have  shown  that  semidiurnal 
internal  tides  are  greatly  enhanced  by  spring  (barotropic) 
tides  and  can  generate  high  frequency  oscillations  with 
marked  peaks  at  periods  of  T=71  and  248IN.  These  deep-water 
internal  tides  could  be  generated  by  the  incidence  of  the 
barotropic  tide  on  the  continental  shelf,  to  -hen  propagate 
offshore  normal  to  the  coast  line  [Rattray  1950].  Their 
amplitudes   would   then  be    given    by: 


where 


tf-      ^ 


and    "A"    is   -.he    amplitude    of      the    internal    wav-,       "2H'1    is   the 
tidal   range,    !,h"    is    the   iepth   to    -he    thermocline,    "D"    is   the 

deep-war er    depth,         "d"    is    the      depth    of   -he      shelf     (assumed 


constant) ,    "W"    is   the    width    of    the    continental    shelf, 


IIt-H       •} 


the    Coriolis 


the  frequency  of  the  internal  wave,  "f"  is 
parameter,  njn  is  the  acceleration  due  to  gravity,  and 
"  ix.  "  is  the  relative  density  difference  between  the  two 
layers.  Hecht's  [1971]  observed  amplitudes  exceed  those 
estimated  by  the  above  formula  (23H  obsevei  vs.  15. 3M 
computed  at  spring  tide).  On  the  continental  shelf  (U6°N 
003°50W),  tearbottom  current  recoris  occasionally  exhibit 
bursts  of  periodic  oscillations  (ca.  30MIN)  followed  during 
a  few  hours  by  apparent  white  noise  fluctuations  [ Le  Floch 
1969  ]. 
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17.     GEOPHYSICAL    SITUATION    OS    AUGUST    20    1978 

A.        HETE0R0L03Y     AND    WIND    STRESS 

On  the  23th  of  August  1978,  a  high  pressure  ridge 
extended   from     the    Azores   to      Western    Europe.  The   average 

surface  pressure  over  the  Bay  of  Biscay  was  1023MB  (Fig. 
4.  1)  .  A  strong  stcrra  generated  on  the  15th  of  August 
offshore  Labrador  was  dying  over  Iceland  (Fig.  '4.3).  A  weak 
cold  front  was  3  00NM  west  of  the  Bay  of  Bisoay  or.  the  19th 
at  1200  GMT  and  probably  underwent  frontolysis  in  the  area 
of    interest. 

From  ship  reports,  the  evolution  of  the  in-situ  atmos- 
pheric situation  (ca.  470N,005°W|  can  be  analysed  (Table 
IV).  Between  200300Z  and  200600Z  the  cold  front  was  present 
near    the    location   of    observations    ao    45O40N,    304°30W. 

TABLE    IV 
METEOROLOGICAL    0BSSRVATI3SS    ca.47N,005W 


I    Day-Hour 
I           GMT 

Surface 

T .           ~-              T- 

|  Wind    from | 

Dry    =iz 

—I" 

I 

Moist 

a  ■ 

— r- 

Clouds 

pressur  = 

\  Speei    M/2 

i 

tea?. oc 

I 

i 

temp. 

°c" 

! 

_  i 

oct  an  ts 

I     191200 

1018.2 

I    250/5 

I 

21.3 

T 

! 

18.  5 

\ 

a 

I     191500 

I 

I 

)     191800 

1 

I 

|     192100 

i     1018.0 

I    25  0/4.5 

1 

20.2 

I 

18.  3 

5 

|    200000 

1020.0 

I     200/3 

I 

18.5 

! 

17.5 

0 

!    200300 

1019.0 

j    270/2  .5 

1 

17. '4 

! 

17.0 

7 

|    200600 

1020.3 

I    000/5 

1 

18.5 

I 

17.  3 

8 

!    200900 

1 

! 

I    201200 

1022.3 

|    030/6 

I 

19.5 

! 

17.2 

6 

I    201500 

1023.0 

!    000/2.5 

1 

18.7 

I 

16.  1 

5 

|    201800 

1022.0 

|    000/6 

| 

18.2 

I 

16.  0 

2 

I    202100 

1021  .0 

!     010/7.5 

I 

17.5 

I 

16.  1 

3 

I    21000  0       I 

1321  .3 

I    010/7 

1 

17.5 

I 

15.  4 

4 

I     21030  0 

\ 

I 

I    210600 

1323.3 

I     040/5 

\ 

17.5 

t 

17.5 

i 

2 

L 
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From  syrootic  charts  and  these  ship  rspor-s  a  comparison 
can  be  made  to  select  the  SEASAT  sea tterometer  winds  a*  the 
time  of  the  SAR  iaage  (fig.  a. 2).  lit  2Q2142Z,  a  local  high 
pressure  center  was  in  the  vicinity  of  480N,  010°  w,  '-has 
establishing  an  anticyclonic  surfacs  flow  decreasing  from 
ca.  10M/S  ssuth  of  Penmarc'h  (north  of  the  SAR  image)  to 
ca.  5  M/S  a*  45°N  006QW  (south  of  tis  SAR  image)  .  The  best 
estimate  for  the  direction  is  a  flow  towards  200°T  (winds 
north-northeast)  . 


B. 


SWELL 


A  heavy  (mean  significant  wave  height  of  2.5M),  and  long 
(mean  wavelength  of  230M)  swell  had  originated  2000KM  away 
(Fig.  4.U)  and  was  currently  prevailing  over  the  sea  state. 
It  was  propagating  roughly  130OT  (froa  shipboard  meteorolo- 
gical reports)  .  The  generation  areas  were  estimated  by 
atmospheric  charts  (Fig.  '4.3)  tc  bs  in  the  southwest  quad- 
rant of  the  surface  pressure  low.  Distances  and  time  of 
propagation  with  respect  to  the  SAR  inage  location  and  time 
yielded  estimates  of  the  wavelengths  that  could  be  expected 
in  the  Bay  of  3iscayr  August  2  0f  as  a  function  of  m ear- 
direction.  The  deep  water  surfacs  gravity  waves  equations 
used    are   the   following: 


the   period   is   expressed    as 
the   group   speed    is  given    by 
thus    yielding 


t=\/!zE 

9 


c,-  «  v/3E 

J  2      ^    2-TT 


ir 


9        t2- 

with    "D"  ,      distance    from    t h e   generation      area    and    "fc",      time 

elapsed   since  the   oriainal    wind    generation. 

No  swell  spectra  were  available  to  compare  with  the 
patterns  visible  :n  the  S&R  imager  but  a  correlation  with 
significant  wave  heights  as  measured  by  the  altimeter  of 
SEASAT  was  possible.  The    most   interesting    feature      in   -he 
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altimeter  pass  of  210721Z  (ascending  node  passing  by  U6°N, 
005°w)  was  a  sudden  increase  in  significant  wave  height  from 
an  average  of  2M  to  mors  than  4  1  over  a  short  distance 
(15KM)  at  position  47015'y,  006°00'tf  (Fig.  4.5),  i.e.  near 
the  shelfbreak.  This  astonishingly  abrupt  variation  could 
hardly  be  related  to  any  atmospheric  mesoscale  process.  It 
will  be  shown  (Cf  VI-E)  that  it  may  be  the  consequence  of  a 
strong  interaction  with  current  shear  occurring  at  the 
shelfbreak.  According  to  sailers  experienced  with  this  area 
[Chiari  1982],  crossing  the  200H  isobath  is  often  accompa- 
nied by  a  characteristic  ohange  in  sea-state.  An  historical 
reference  is  given  by  Le  Dancis  [1943],  who  attributes  the 
denomination  of  La  Chapelle  Bank  to  an  observation  mad£  by 
the  "Sieur  da  La  Chapelle"  in  1695  of  surface  turbulence 
wrongly  associated  with  extreme  shoaling.  This  apparent 
shoal  was  resighted  in  1754  and  a  rock  was  marked  on  the 
charts    for   the    following    decades.  A    well    known    example    of 

such  a  current-swell  interaction  is  the  occurrence  of 
"giant"  wav=s  [Huthnar.ce  1980]  near  the  shelf  edge  off  the 
East  coast  of  South  Africa.  These  are  thought  to  arise  at  a 
caustic  associated  wi^h  the  refraotion  of  surface  waves  by 
the  Agulhas  Current.  In  the  same  way,  an  alongshore  current 
at  the  shelfbreak  in  the  Celtic  shelf  could  explain  the 
increase    in    wave   height    observed   there. 

C.       MEAN    AND   TIDAL    CURRENTS 

Data  available  for  iirec":  estimation  of  the  mean  curren^ 
and  its  variations  over  the  area,  are  proviiei  by  the  ieeo 
water  subsurface  flooring  at  U6°57'tf,  009°57MW  with  current 
meters  a1:  1 52  U ,  3032  and  3976a  (350O  off  the  5AR  ground 
xrack  and  180K*  offshore  of  the  shelfbreak)  (Cf  1-E).  Lew 
fraquency  currents  at  the  three  levels  do  not  exhibi*  any 
obvious      correlation    and      vary   over      all      directions    with      a 
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time-scale   cf  three    months.         The      n  =  an    currants    in    the   last 
10    days    before    the    20th   of    August    ace    r aspect ively: 

at    1524M    toward   010°T    at    0.05M/S 

at    3032M    toward   27Qor    at    0.034H/S 

(veering    33D°T   in   last   3    days) 

§.1  19  76M    toward   315»r    at    0.001M/S 
On   a    shorter  time   scale,      progressiva    vector    diagrams    show   a 
semidiurnal    pariodicity      at    each    of      the    three      levels     (Fig. 
4.6)  ;    the   mean    speeds  are    : 

§.*    lUll*    G.08M/S 

at    30  32M    0.05M/S 

at    39  76M    0.00  15M/S 
At    20    21    40,      time    of  tha   SEASAT    passage,      the   instantaneous 
current    vectors    were    : 

§.1  1524M    toward   315°r    a4!    0.08.1/S 

at    30  32M    toward    18 5° r    at    0.06.1/S 

al  3976M  toward  270°!  at  0.0Q15H/S 
The  passage  of  SEASAT  on  20  August  followed  a  spring  tide 
ranging  110  on  the  Frenoh  coefficient  scale  [SHOM  1968] 
(i.e.,  amplitudes  1 535  greater  than  tha  average  spring  tides 
of  coefficient  95).  At  the  time  of  tha  imag-,  the  surface 
tide  at  the  shelfbreak  was  at  high  tid=  plus  four  ar.d  a  half 
hours       (Table      7).  The   surface      currents      were      estimated 

towards    230    or    at   0.85M/S    (Table    4.5),    [SHOM    1968]. 

D.        SPATIAL    THERMAL    STRUCTURE 

The  thermal  structure  over  the  irea  was  established  from 
a  set  of  69  X31V5  (Cf  1-8).  The  time  span  of  this  tempera- 
ture sampling  was  from  09  00  17  August  to  0500  24  Auoust. 
The  spatial  Limits  are  between  43°  aid  '48° North,  and  10° 
West  and  the  coastline.  The  distribution  of  this  sampling 
is  conveniently  centered  around  tha  radar  image,  in  partic- 
ular   in      the   vicinitv      of    the    continental      slope    and      with   a 
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TABLE    V 
CDEFFICIENTS    MID    TIHES    DF    HIGH    TIDES 


Days    af 

In    BREST  | 

At    Shelf] 

Coef  f 

j  Amplitude  I  Amplitude! 

Auaust 

GMT    Hour| 

break         | 

(in    Brest:    |  at    Break    l 

15 

00    33 

00    13       | 

60 

|     1  .  9  3  5          |     1  .  2  5  M 

13    09    ( 

12    49       | 

D  3 

I     2.12             I     1.38            I 

16            | 

0  1     42     J 

01    22       I 

73 

|    2.34            |     1.52            | 

14     12 

13    52 

33 

I     2.57            |     1.67            | 

17 

0  2    40     | 

02    23       { 

37 

|     2.79            |     1  .81            I 

15    06     I 

14    4  6       | 

94 

|     3.04             |     1.97            | 

18            I 

33     32     I 

03    12       | 

100 

I     3.20            I     2.03            I 

15    56    | 

15    36       | 

134 

!     3.33            I     2.16            | 

19 

34     19     | 

03    59 

133 

I     3.46            J     2.25            | 

16    42    I 

16    22       | 

113 

I     3.50            |     2.28            I 
,.  i _. j                         j 

20 

35    05     | 

04    45        | 

113 

j     3.50            j     2.2  3            j 

17     26     | 

17    06       | 

139 

|     3.48            I     2.26            | 

SAR    imag* 

■:    at    2  142Z 

and    at    the 

sheLfb 

reak    High    Tide   +  4.5H    I 
_i _.  _ i ..           i 

21 

35    47 

05    27 

13S 

i'                       i                       i 
!     3.41            |     2.21            | 

13    03     | 

17    43        | 

132 

j     3.26             |2.11            I 

22           | 

06     28    | 

06    03       | 

97 

j    3.12            I     2.02           | 

18     49     | 

18    29        | 

93 

I     2.90            !     1.83            I 

23           ! 

0  7    09     | 

06    4  9       I 

33 

2.60            |     1.69            I 

I 

L I 

19    30 

19    10 

76 

|     2.45            |     1.59            | 

maximum  density  on  the  eastern  half  of  the  ground  coverage, 
as  shown  by  a  plot  of  the  spatial-t emporal  distribution 
(Fig.  4.7)  with  respect  to  the  distance  from  the  shelfbrsak 
(y-coordinate  defined  in    1-C   and    Fig.    1.1). 

1  •      Analysis   of    the    Mean 

A  smoothed  transect  was  established  by  averaging 
along  the  x-axis  with  a  "  y"  step  increment  of  20O  (Fig. 
4.3)  .  The  isotherms  generally  sloped  down  from  the  shelf 
region  to  the  open  sea  (on  the  average  10M  over  100KM.).  A 
marked  mixir.  a  was  apparent  in  the  region  of  the  shelf  break 
(-10<Y<10)  with  shoaling  :f  the  isothermal  layers  above  the 
thermociine  (14°C)  and  sinking  belaw  it.  Dynamic  heights 
relative  to  the  smoothed  transect  yield  a  mean  surface  gecs- 
^roohic  flow  relative  to  150M  of  the  order  of  .2M/S  and 
about    20KH    vide,    towards   the   northwest   at    the    shelfbreak. 
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TABLE    71 
SURFACE    TIDAL   CURRENTS    AT    * 7 o 2 «* •  N    006005'? 


!       Hour 


Direction | 


GMT 


ot 


Spr  ing 


Speeds 
Neap 


/s 

T 


AQGUST     20 


f 1 

OH    45 

—     ■  ■-—■ —  \               — 

,  i I  .High   Tide 

I       05    05 

120 

0.  25 

I            0.15                      0. 3 

!       06    05 

190 

0.  4 

!            0.2               !          0.45 

I       07    05 

20  0 

0.  5 

I            0.25            |          0.55 

I       OS    05 

210 

0.  6 

I            0.3               |          0.65 

|       09    05 

22  5 

0.  75 

I             0.4               |          0.85 

1       10    05 

23  8 

0.  75 

I             0.4               |          0.95 

!       1  1    05 

23  6 

0.  5 

I             0.25             |          0.55 

I       1126 

25  0 

0.  4 

I            0.2              I          0.U5 

I       12    26 

23  5 

0.  2 

(            0.1               j          0.2 

I       13    26 

10 

0.  3 

I             0.15             |           0.3 

!        14    26 

35 

0.5 

I             0.25             1          0.55 

I       15    26 

50 

0.  5 

I             0.25            1          0.55 

|       16    26 

70 

0.  45 

1             0.25             1          0.5 

I        17    06 

H^ah   T^de 

|       17    26 

120               I 

0.  25 

i            0.15                      0.'3 

I       18    26 

190 

0.  4 

|             0.2               |           0.45 

<       19    26 

20  0 

0.  5 

I            0.25            I          0.55 

I       20    26 

210 

0.  6 

I            0.3               I          0.65 

1       2  1     26 

22  5 

0.75 

0.4                          0.85 

|SAR    image 

—  at   214. 

>   —    Awards 

2303    at    0.85    S/S    

I       22    26 

23  8               | 

0.  75 

I             0.4               |          0. 95 

I       23    26       ! 

236 
25  0 

0.  5 

I            0.25            I          0.55 

I        23    47 
i _.  ._    i 

3.  4 

|             0.2                |           0.4 

The  geometric  approximation,  of  the  shelf  by  a 
s-raighr  line  in  the  averaging  proo=3s,  and  probable  varia- 
tions in  time  and  space  of  this  mixing  region,  reduce  the 
more  dramatic  features  that  can  bs  se  =  n  in  Darticular  non- 
averaged  transects.  The  mean  vertical  -smperature  profile 
is  characterized  by  a  mixed  layer  at  13°C  limited  by  a  sharp 
thermocline  (on  the  average  0.2°C/M>.  On  the  shelf,  a  bottom 
mixed    layer   of    12°C    was    present. 

Taking  an  arbitrary  "level  of  no  motion"  at  150M, 
dynamic  heights  were  compared  for  Bach  X3T  soundings  and  a 
horizontal  smoothed  contouring  was  Bade  (Fig.  4.14).  The 
surface  gecstrophio  speeds  corresponding  to  this  analysis 
are  cyclonic  around  a  stretched  aiiipse  with  a  major  axis 
coinciding  with  the  continental  s1od=.  Nor-.ti*est  currents 
on      the    uooer      Dart      of      the   slooe      ars      of      .5&/S    at      ~heir 
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maximum,      vs.      a   returning    southeast      flow   of    ,23/S    over   ths 
ri  s  9 . 

A  vertical  modal  analysis  (Cf.  VI-B)  ,  yields  baro- 
clinic  internal  Rossby  radii  of  deformation  of  4.5KM  on  the 
shelf   and    25KS    in   the  deep    ocean.       (Cf    VI-B) 

2.  An  a  1  ys  i  s    of   the   Variability 

The  spatial  variability  of  the  thermocline  depth  was 
computed  using  a  subset  of  the  total  X3T  data,  which  was 
selected  to  avoid  unrepresentative  (spares  and  dispersed) 
values.  The  variance  of  the  depth  of  a  given  isotherm  was 
calculated  in  each  rectangular  bin  Df  10:<M  width.  Standard 
deviations  of  the  1 3°C  isotherm  depth  from  the  above  means 
were  computed  (Fig.  4.9)  .  The  strong  increase  in  the  values 
at  the  position  of  the  idealized  shelf  break  could  either 
mean  a  significant  increase  in  the  time  variability 
(internal  waves),  or  express  the  fart  that  the  mixing  is 
taking  place  over  a  very  short  distance  compared  with  tha 
width    (10KK)    of    the    bins    used    for    averaging. 

3.  Select  ad    Reorese  nta  tiye    Transects 

Seme  specific  transects  were  then  considered  to  give 
a  quasi- instantaneous  spatial  image  of  this  mixing  process. 
One  interesting  subset  of  BT's  was  acquired  by  a  transiting 
ship  sailing  towards  Ouessar. t  (Ushant  Island)  and  whose 
track  crossed  the  shelfbreak  at  47033N,  D6<>H  on  August  23.  A 
sampling  rate  of  one  per  hour  (i.e.  30KM  spacing)  enables 
drawing  a  quasi- synoptic  transect  parallel  to  the  SAR  swath 
but  60KM  west  of  it  (Fig.  4.10).  Other  transects  were  made 
by  selecting  3T's  close  in  time  ani  conveniently  spaced  in 
the  vicinity  :>  f  a  perpendicular  to  the  isobaths.  A  se~  of 
three  transects  is  depicted  on  a  geographic  grii  (Fig.  4.11) 
with  other  non-synoptic  ^ransects  reproduced  from  various 
sources     (Fia.      4.12).        The    scaling    adopted    is    identical    for 
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all  and  bas^d  on  the  "alon g-^rack"  distance  from  the  shelf- 
break.  Ccmmon  characteristics  appaar  to  be  a  narrow  pelagic 
front  anchored  on  the  sea-floor  at  the  upper  part  of  the 
shelfbreak  (ca .  180.1  depth}  ,  with  01=  or  two  upwelling  zones 
adjacent  to  this  narrow  intermediate  mixed  layer.  The  ther- 
mocline  is  stronger  .on  the  shelf  than  in  deep  waters  (Fig. 
4.13).  In  chapter  VII,  an  idealized  3-layered  ocean  was 
extrapolated  from  these  -elected  transects  and  the  corre- 
sponding geostrophic  currents  were  compared  to  the  observed 
circulation. 

E.       TEMPORAL   THERMAL    ANALYSIS 

Daily  averages  of  the  isothermal  depths  (Fig.  4.15)  show 
a  significant  heating  in  the  upper  layers  on  20  Auaust  and  a 
general  broadening  of  the  thermocline  between  17  and  24 
August.  Assuming  a  propagation  offshore  from  the  shelf  of 
tidally  generated  waves,  diagrams  have  b=en  established 
displaying  foe  a  given  isotherm,  variations  of  its  depth  as 
a  function  of  time  and  space  (Fig.  4.16).  These  variations 
in  depth  were  obtained  by  subtracting  the  mean  along  the 
x-directicn  to  eliminate  any  permanent  trend  and  enable  a 
contouring   of      relative    values.  Da?    spatial      axis    is      the 

along-track  (y-coor dinat e)  distance  from  the  shelfbreak;  ^he 
time  axis  in  hours  has  bsen  "folded  ov=>r"  two  semidiurnal 
(12. 4H)  periods.  Henc=  all  tins  coordinates  are  *"aken 
within  24  hours  by  subtract  ion  (or  addition)  to  the  actual 
time  of  observation  of  an  adequate  nultiple  of  the  semidi- 
urnal   period. 

No  contouring  algorithm  could  ba  applied  to  such  sparse 
and  greatly  contaminated  data  (by  various  phenomena:  vari- 
able wind  stress,  radiation  heat  exchange,  even  high 
frequency  nonlinear  internal  waves).  A  hand-contouring  was 
attempted  under    the    fcllowincr      oriicloles:       one    contour   line 
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was  draws  fcr  a  z?ro  relative  depth,  another  line  followed 
the  maxima  and  the  last  one  followed  the  minima.  This  was 
done  for  three  temperatures  (12,14 ,  15°C)  .  Each  of  '-hose 
diagrams  exhibits  a  constant  pattern  obliquely  oriented  on 
the  axis,  and  thus  corresponding  to  a  propagation  of  the 
variations.  The  slope  of  those  patterns  (distance  of  propa- 
gation during  a  given  interval  of  time)  was  found  to  be 
approximately  1 . 1M/S  (equivalent  to  0.9H/S  perpendicular  to 
the  shelfbreak)  .  A  modal  analysis  described  in  VI-B  showed 
that  high  frequency  internal  waves  of  period  90MIN  would 
have  this  same  wave  speed.  Furthermore  the  phase  of  -his 
propagation  implied  a  maximum  deepening  event  a*  my  =  0"  a1: 
about  1900Z,  20  Add.  Conversely,  at  the  time  of  the  SAR 
image  the  crests  of  these  propagating  oscillations  were 
found  to  be  it  Y=10  and  Y=60.  Maximum  values  of  relative 
depth   excursions    show  general    increasing    values    with   dep-h: 

a.  Fcr   isotherm    12°C    mean    dep^h    of    35M    with    55M    excursions 

b.  For   isotherm    14<>c   mean    depth   of   45M   with    15M    excursions 

c.  For  isotherm     160C    mean    depth   of    30M    with    1  7M    excursions 

At   deeper     depths,      the      NEADS    mooring      temperature   data 
allows      estimation      of      the    vertical      isotherm      movement      by 
relating      semidiurnal  temperature      fluctuations      at    a      given 
depth    to    the   ol imatolcgical    temperature    profiles.       Data    from 
the    current   meter  pressure    gauge    showed    that    vertical   oscil- 
lations   of   ^he    array    itself    were    negligible. 
At    1500M,      teiperature  variations    of    ,3°c    in   a   mean    gradient 
of    3    10"  3    oc/M    yield    vertical    amplitudes    of    about    20M. 
At    3000M,    temperature   variations    of    .  15°C    in    a   mean    gradient 
of    2    10~*    oc/H    yield    vertical    amplitudes    of    more    than    U 0M . 
At      3500M,        temperature      variations    are      at      the      limit      of 
discrimination      of    the     thermal      gaige      but    reveal      vertical 
movements   of  at    least  20M. 


58 


F.        SUMMARY 

SEASAT  flaw  over  "rh0  shelfbreak  3  f  the  Bay  of  Biscay  a4: 
2142Z  (after  sunset)  on  20  August  1973.  A  high  pressure  cell 
had  dominated  tha  area  since  early  in  the  morning,  bringing 
a  southward  wind  of  6tt/S.  Tha  cloud  cover  gradually 
decreased,    and    the    air  temperature    da-creased    to    17oc. 

A  well  established  swall  was  entering  tha  Bay  from  the 
northwest,  nearly  tangential  to  tha  isobaths  in  the  vicinity 
of  the  shelfbreak  where  a  significant  increase  in  wave 
heights    was    observed    from   the    SEASAT    altimeter. 

The  tide  was  a  spring  tide  at  high  level  plus  5  hours  a*. 
the    shelf break.  The   ebbinc    surfaoa      current    on      the    shelf 

side    was    at    its    maximum    toward    2  4  0°  r    at    0.8  5:1/5. 

The  ocean  was  strongly  stratified  with  sea  surface  Temp- 
eratures of  19°C,  a  thermocline  gradient  of  0.3°C/M  at  an 
average  dep^h  of  30M,  and  a  bo^ooa  tidally  mixed  layer  on 
the  shelf  of  12°C.  The  first  moda  internal  Rossby  radii  of 
deformation  were  of  4.5K3  on  the  shelf  and  of  25KM  in  deep 
wa  t  e  r . 

A  spatial  thermal  analysis  showed  a  marked  fron4:  a4:  the 
shelf  break  with  intense  mixing  and  upwelling.  There  was  a 
westward    surface   geostrophic   jet    at    the    snelfbreak. 

A  temporal  thermal  analysis  gave  some  evidence  of  an 
internal  semi-diurnal  "disturbance"  propagating  offshore 
from  the  break  with  a  speed  of  0.9S/S  and  with  the  last 
crest  in  the  vicinity  of  "Y=10"  it  the  time  of  the  SAR 
image.  This  propagation  sp^ed  being  significantly  smaller 
than  tha  computed  phase  speed  of  the  internal  first  modes 
tides  (3.8  and  1.7M/S),  either  some  higher  modes  are  presen*: 
or  high  frequency  waves  are  generaoai  by  the  tide  but  even- 
tually decoupled  from  it  and  trapped  in  the  vicinity  of  tha 
thermocline  (where  the  largest  Brunt- Vaisala  frequencies  can 
be    found)  . 
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▼•    SAR    IMAGE    ASALJSIS 

A.       ORBIT    PARAMETERS    AND    GEOGRAPHIC    REFERENCES 

The  auxiliary  grid  introduced  previously  (1,4)  will  be 
used  to  cor. valient ly  locate  features  with  respect  to  the  SAR 
image  and  tc  tha  average  shelfbreak  contour.  The  geographic 
fit  of  the  SAR  image  is  estimated  ~d  be  accurate  to  a"-  least 
+/-    5KM   in    the    vicinity    of    the   shelf-break. 

At  2  142Z  Aug  20,  SEAS  AT  was  3  a  a  descending  r.oda  over 
Franca  with  the  SAR  ant  anna  oriented  to  the  northwest. 
Consequently,         a    positive       (northward)  shear    of      eastward 

current  would  appear  on  the  image  as  artificially  compressed 
with  enhanced  Bragg-scatter  (Cf  Fig.  2.3)  In  other  words,  a 
current  directed  toward  the  northwest  would  appear  on  the 
SAR  image  with  a  bright  northeastern  boundary  and  a  dark 
southwestern  one. 

TABLE    711 
SEASAT   ORBIT    785    PARAMETERS 


| Ascending   node    at 
I  East    Longitude 

I  Descending    node    at 
I  East    Longitude 

I  Argument    of  Perigee 
[Eccentricity 


Orbit    785    on    August    20    1973    (Julian    Day    232) 


2103. 55Z 

168.43    32     (decimal    dearies) 

2154. 1 4  Z 

335.39    58     (decimal    daqrees) 

w=140 

e=o.:ms 


j 
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B.        HEURISTIC    ANALYSIS    OP    THE    SAR    IHAGE 

The  patterns  visible  in  tha  optically  processed  SAR 
image  may  be  grouped  in  four  distinct  categories  according 
to  their  geographic  location  and  to  their  geoaeoric  charac- 
teristics. They  are  discissed  separately  under  the  titles: 
"inner  shelf",  "outer  shelf",  "generation  area",  and  "propa- 
gation area".  The  image  discussed  being  a  positive  print, 
the  lowest  baokscatter  radar  respois?  will  be  referred  to  as 
dark  areas  and,  conversely,  the  enhanced  signal  returns  will 
be   referred   to    as   bright   or    light   ir=as. 

A  correlation  between  image  characteristics  and  ba thy  me- 
tric features  is  made  and  must  taks  into  account  two  inaccu- 
racies. First  the  error  in  mapping  the  SAR  image  is 
estimated  as  5KM  at  most.  SeconiLy  the  bathymetric  chart 
used  here  (Fig.  5.5)  was  constructed  by  interpolating  sound- 
ings of  different  sources,  with  gaps  of  ca .  10KM  in  some 
places  of  the  continental  slope.  Hance,  the  relative  posi- 
tions and  shapes  of  the  canyons  drawn  are  subject  to 
caution  within  the  same  order  of  magnitude  as  the  SAR 
feature    locations.  Nevertheless,       in    the   case      of   several 

assumed  correlations  between  well-defined  imaae  and  bathy me- 
tric features,  the  confidence  nay  be  increased  by  the 
perfect   relative    fit    of   each   of  the    points     (Fig.    5.4). 

1 .      Inner  Shelf 


(Y>40)  -    In      this    area,         mai  i  fe  stations    of 


n  t  e~  na  1 


activity  are  of  low  contrast  and  Lark  any  obvious  organisa- 
tion. Separated  curved  lines  can  nevertheless  be  identified 
with  an  average  radius  of  curvaturs  of  12  K?!.  Patterns  of 
similar  shapes  can  be  associated  at  spacing-  of  280  (Fig. 
5.2) ,  but  thair  relation  to  any  kind  of  periodic  phenomena 
is  not  proven.  They  could  be  due  to  mesoscale  atmospheric 
variability,  or  linear  slicks  associated  with  a  localized 
convergence    (Z  f    oara    II-C-2)  . 


77 


The  lowest  limits  of  wind  speeds  suff  icier.4-  for  a 
radar  return  are  attained  in  the  northern  par-  of  -his  area, 
particularly  south  of  the  "Chausses  of  Sein"  and  leeward  of 
the  Britanny  3oast  where  mesoscale  land-sea  breeze  may  well 
have  been  predominant.  (A  ship  saiLing  100°r  can  clearly  b^ 
seen  with  its  wake  .  virtually  displaced  according  to  the 
Doppler    effect    described    in    II-B-1). 

2-      Outer  Shelf 

(10<Y<40)-  Internal  wave  group  patterns  are  apparent 
in  this  region.  They  appear  as  arcs  of  circles,  concav<=  tc 
the  west  with  a  variable  radius  of  curvature  of  minimum 
value  15KM.  They  are  United  to  the  north  by  a  caustic  in 
which  they  merge;  this  line  is  brigh'-Br,  i.e.,  of  greater 
backscatter.  They  are  broken  fc  the  south  by  another  line 
of   brighter   characteristics    still.  These   patterns   seem   to 

exist  over  the  entire  width  of  tie  3AR  imacre,  bu4:  "hey 
appear  with  more  contrast  on  the  western  par-  of  it.  Som= 
periodicity,  of  about  29KM  length  op.  the  x-axis,  car.  be 
detected . 

The  wavelength  of  these  patterns  is  variable  in  two 
directions.  (1) toward  the  so-callad  caustic  to  the  north, 
the  wavelengths  decrease  to  zero,  starting  at  a  maximum 
value  on  the  southern  boundary.  (2> toward  the  southeast 
(increasing  x-direction)  ,  which  is  thair  apparent  direction 
of  propagation,  the  average  wavelengths,  measured  at  mid- 
distance  between  the  two  boundaries,  is  of  the  order  of  1K!1 
on  the  western  side,  increasing  to  2KJ1  in  the  middle  of  the 
swath. 

The  most  conspicuous  waverrests  are  located  at 
x=-25,  where  a  high  level  contrast  particularly  enhances 
those      curved   lines.  Each    of      thsse    arcs      appears    as      two 

narrow  contiguous  stripes:  a  iark  (low  backscatter)  one  on 
the    concavity      and    a    bright       (high    backscatter)       one      on   the 
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convexity.  Hence,      over      one      wavelength      on    a      direction 

perpendicular  to  the  crests  and  following  tha  assumed  dirsc- 
tion  of  propagation  (from  the  inferior  to  the  exterior  of 
the  arcs),  the  backscatter  varies  approximately  in  tha 
following  way:  73^  of  the  length  at  a  uniform  level,  2 0%  at 
low  backscatter  level  and  1 0%  at  a  iigh  level  on  the  leading 
edge.  According  to  the  hydrodynamical  theory  of  SAR  imaging, 
this  result  means  that  tha  internal  waves  responsible  for 
these  pattens  are  long-crested  with  deep  narrow  troughs 
(Fig.    2.  4)  . 

The  southern  boundary  is  a  general  line  brighter 
than  the  nortaern  caustic  particularly  on  the  wes-ern  part 
of  the  image,  where  it  is  conspicuous  both  by  its  brighter 
characteristics  and  by  tha  obvious  discontinuity  it  creates 
in  the  general  patterns.  This  bouiiary  is  located  approxi- 
mately (♦/-  5KM)  over  the  180M  isobath,  thus  corresponding 
to    the   upper  break    in  the    continental    slope. 

3.      Continental    Slope   or.  [ianeri^ion    Area 

(0<Y<10)-  This  area  covers  the  upper  part  of  the 
slope  between  tha  150  and  1000M  isobaths.  Once  again  tha 
mesoscale  features  are  concentratad  on  the  westarn  side  of 
tha  image.  Tha  average  backscatter  is  weaker  (darker 
shading)  than  on  both  adjacent  areas  (shelf  and  deep  ocean). 
The  northern  limit  is  tha  band  of  lighter  characteristics 
described  above. 

The  limit  to  the  south  is  (conversely)  a  discontin- 
uous band  of  /ery  low  backscatter  about  1 5 0 0 M  wide,  which  is 
clearly  visible  from  x=-60  to  x=-23  aid  from  x=-5  to  x=15. 
The  straightnass  of  this  line  (vs.  tha  other  patterns)  and 
tha  fact  that  it  is  not  bordered  bv  any  alternating  clearer 
stripes,  seem  to  indicate  that  it  is  due,  not  to  a  propa- 
gating phenomenon,  but  to  a  stationary  (if  not  steady) 
internal    feature.  Thasa    could      ba    symDtoms      of    a      surface 
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divergence      according     to      the      kinematical      theory      of      S&H 
imaging. 

In  this  area,  trell-organi ssd  coherent  circular 
patterns  suggest  cylindrical  spreading  of  internal  waves 
from      several  distinct      point    sour:5S.  Assuming   -hat      the 

direction  of  propagation  of  a  givei  wavepacket  can  be  esti- 
mated both  by  the  curvature  of  the  :ras:s  and  by  the  rela- 
tive wavelength  dispersion  sense,  the  following  main  point 
souces  can  be  identified  ind  given  names  for  further  re: 
e  n  c  e  . 


<=+■  ar  - 


P1    centers!    a"         X  =  -42 

Y=3 


=    1*70(3  ON 

=D05O30W 


P2    centered    at        X=-25 

YO 


=    46052N 
=D05°18W 


?3    centered    at        X  =  0 

Y  =  5 


=    U6045N 
=005<502W 


?U    centered    a4*         X=4  5  =   46  0  25N 

Y  =  3  =004°42W 

Point  P1  is  at  the  head  cf  Penmarc'h  Canyon  at  300M 
depth. 

Point  P2  is  at  the  head  cf  a  narrow  and  oblique  but 
deep,    unnamed  canyon. 

Point  P3  is  at  the  head  of  3uilvinec  Canyon  at  BOOM 
depth. 

Point  P4  is  at  the  head  of  3lavet  Canyon  at  300M 
depth. 

It  is  observed  that,  at  the  Oiet  canyon,  there  is  no 
evident  corresponding  internal  wave  pattern;  a  plausible 
reason  for  this  could  be  that  this  oinyon  thoujti  steep,  does 
net    go    up   to  the   shelfbreak    itself    (Fig.    5.5)  . 
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At  ?1  the  leading  wave  cf  the  nearest  group  is  15^^ 
from  the  center.  Wavelengths  from  the  leading  trough 
decrease  from  2KM  to  0.7KM  with  a  range  of  +/-  50*.  Som  = 
interference  02.11  be  seen  with  a  second  circular  pattern  of 
poorer   definition  centered    around   ?'1    (X=-44, y=- 1 0)  . 

At  P2,  two  wave  groups  propagate  in  two  distinct 
directions.  The  first  group  propaga-es  toward  130°?  along 
the  X-axis  with  its  leading  wave  it  10KM  from  the  center. 
The  second  group  propagatss  with  a  33°  inclination  offshore 
(160°T),  its  leading  wava  being  at  18KM  from  the  center  and 
having   the   greater    wavelength   in    tha    group    as    above. 

A*:  P3,  no  obvioas  group  ran  be  identified  bu*-  the 
pattern  appears  as  a  superposition  of  faint  Ellipses  with 
their    major   axis   oriented    130°T   along    the   X-axis. 

In  the  vicinity  of  P4,  patterns  ara  still  mora 
confused,  but  a  qaasi-cirr ular  group  axtends  around  it  at  a 
distance    cf    17KM. 

U.      Continental    Riss    or    Prooaaat i oa    Area 

(-9  0<Y<0)  The  surface  pattarns  here  ara  more  blurred 
than  in  the  three  areas  ^ascribed  above.  Features  are  on  a 
larger  scale  and  do  not  present  the  sane  east -west  vari- 
ability. The  brightness  i=  uniform  over  the  whole  area.  The 
general  geonretric  organization  is  that  of  arcs  of  circles 
roughly  centered  on  the  point  sourras  lefined  above,  or  more 
generally,  on  the  1000M  isobath.  Some  alternating  black  and 
white  bands  may  be  identified,  buc  ara  larger  and  do  not 
contrast  sharply  on  a  grey  backdrop  (in  contrast  to  the 
so-called  generation  area).  Tha  distribution  of  those 
circular  bands  is  not  always  clear,  but  some  particular 
patterns  with  similar  shapes  can  ba  found  at  a  spacing  of 
ca.       23    KM. 
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C.  ANALYSIS    OF    THS    SWELL 

The  remotely  generated  swell  is  long  and  well  organized, 
thus  giving  generally  well  identified  patterns  on  the  SAR 
image.  Some  apparent  geographic  variability  in  direction  and 
crest  length  ran  be  observed  and  could  be  indicative  either 
of  variations  in  the  nonlinear  imaging  mechanism,  or  of  real 
deformations  in  the  swell  field  by  interaction  with  surface 
currents  or  mesoscale  winds.  The  wavelength  of  the  swell  is 
237M     (at    U7°N,  004O.  30'  »)  ♦/-    30M.         This    corresponds      to    a 

period  of  12.3  seconds,  a  phase  speed  of  19.2  K/S#  and  a 
group   speed      of    9.6    M/S.  The    direction    of      propagation    is 

well  defined  and  easily  measurable  :n  the  shelf  and  offshor= 
Spain  on  the  southern  par1!  of  the  5\3  image.  In  the  vicinity 
of  the  Armorican  shelfbreak,  the  patterns  are  nore  confused 
and  mix  with  the  internal  wave  field.  On  the  shelf,  the 
general  orientation  is  115°T  (at  47°N, 004QW) .  In  the  open 
ocean  (46<>N ,  006°  W)  ,  it  is  135°T.  South  of  the  shelfbreak 
it  seenis  to  drastically  bend  towards  the  south  as  if 
refracted  by  an  opposing  curren*-.  This  particular  feature 
will    be    discussed    further   in   sectioi    VI-E. 

D.  SOHMAEY 

On  the  shelf,  and,  particularly  in  the  first  30  KM 
inshore  of  -'-he  shelfbreak,  ^he  general  orientation  of  the 
surface  patterns  suggest  an  internal  motion  propagating 
along  the  isobaths  and  toward  the  east.  A  strong  and  local- 
ized divergence  follows  *-he  shelfbreak  at  180M  of  depth, 
clearly  limiting  to  the  north  a  generation  area  on  the  slope 
where  several  point- sour ces  can  be  identified  at  ~he  heads 
of  sharp  canyons.  Internal  waves  are  generated  in  packers 
which  are  clsarly  visible  at  short  distances  from  their 
generation  points.  Wavelengths  of  the  internal  waves,  as 
well    as    the      spacing    between   oackeos,       are      highlv    variable, 


bat  both  increase,  with  the  distance  offshore  from  the  slope, 
The  general  intensity  and  contrast  of  all  patterns  decr=as! 
from  the  northwest  to  the  southeast  (with  decreasing  rang; 
from    the    SAR   antenna). 
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VI.    CORRELATIVE    INTER PRETATI39    OF    THE    SAR    IMAGE 

The  entire  area  of  tha  continental  slope  imaged  by  SAR 
is  dominated  by  surface  manifestations  of  short  internal 
waves.  These  internal  waves  appear  in  groups  which  are  prob- 
ably generated  at  the  shelfbreak  in  association  with  a 
semidiurnal  internal  tide,  but  avaatually  decouple  from  it 
as  both  perturbations  propagate  offshore.  Two  different 
fields      of    internal      waves    have      to      b=    considered;  (1)the 

internal  tidaL  wave  with  a  semidiurnal  period  of  12. 4  hours, 
and  (2)  the  higher  frequency  waves  generated  in  packers  as 
forced  by  tha  internal  tide  interacting  with  shelf  bathy- 
metry, or  by  dissipating  energy  at  higher  wavenumber  by 
breaking  processes.  The  o haracteris tics  of  these  internal 
waves  are  derived  from  their  surface  manifestations  (Cf 
chapter  V).  Zcmparirg  the  results  with  a  synaptic  conven- 
tional analysis  (Cf  chapter  IV)  of  thermal  structure  and 
Fulerian  current  measurements,  gives  a  notion  of  the 
validity  of  an  operational  approaoi  where  a  SAR  imags  could 
be   the   basis   for   hindcasting   -.he    actual    internal    motions. 

A.       TWO    LAYER    MODEL    ON   THE    SHELF 

The  sharpness  of  ^he  thermocline  on  the  continental 
shelf  allows  the  assumption,  in  a  first  step,  of  a  two -layer 
ocean  and  us?  of  the  corresponding  dispersion  relation  to 
compute  phase  speeds  for  given  obsarvad  wavelengths.  The 
dispersion    equation    is: 


*        (W^o+P,^^. 
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Or.  the  ou-'ier  shelf,  for  a  giver,  internal  wave  crest,  the 
observed  increase  of  wavelength  from  the  "caustics"  (Cf 
V-B-2)  to  the  leading  middle  point  averages  25/5.  Working  by 
logarithmic  differentiation  en  the  simplified  form  of  the 
dispersion    equation 


c1  =     AL  L(iiiK) 

J      o 


H 


y 


gives   an    approximate    relation      between    this    wavelength   rela- 


tive   variations    and   th* 


.acteasi 


in    mixed    layer    depth: 


^C  0  U  U  {  k 


U-V\ 


and 


According  to  the  observed  variations  and  to  the  mean 
thermocline  depth  derived  from  XBT  data,  the  thermcclina 
depth  is  estimated  to  vary  on  the  outer  shelf  from  25'1  on 
the    boundaries    to  55M  at   mid-distance    from    each    of   them. 

On  the  coatinental  slope,  a  linear  interfacial  wave  of 
wavelength  1  KM  yields  a  phase  speed  of  0.481/S  whereas  for 
700M  the  phase  speed  decreases  to  3.U7M/S.  Applying  these 
results  to  th?  distance  separating  the  wave  packers  gener- 
ated on  the  shelf,  and  assuming  tie  generation  points  are 
those  bathymetric  features  described  above  (Cf.III-B),  the 
travel  time  would  be  9  hours  for  the  leading  waves,  thus 
implying  a  generation  at  low  tide  with  maximum  ebbing 
surface    currents    (Cf    Table    U.6). 

Conversely,  and  if  the  generation  time  was  known  with 
better  orecision  relative  to  the  tide  phase,  observation  of 
the      position      of      similar         groups      together      with      density 
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gradient      data    could      give    an      estiaate      of   the      therroocline 

depth  . 

This  two-layer  approximation  would  hold  if  the  thermo- 
cline  was  sharp  enough  or  squivalently  if,  at  its  level,  the 
Brunt-Vaisala  frequency  profile  exhibited  a  strong  peak, 
thus  trapping  the  energy  of  high  frequency  internal  waves  in 
an  horizontal  channel.  But  indeed,  the  actual  ocean  has  not 
such  a  simplified  structure  ,  particularly  in  the  "genera- 
tion area",  where  dramatic  mixing  processes  are  taking 
place.  Furthermore,  at  lo»  frequencies  such  as  semidiurnal, 
the  vertical  movements  extend  to  tie  whole  water  column  and 
cannot   be   desoribed    without    a    vertical    modal    analysis. 

B.       NORMAL    MOOES    OF    THE    BAROCLINIC    INTERNAL    TIDE 

Internal  baroclinic  modes  of  the  ocean  were  computed 
from  the  density  profiles  extracted  from  observed  temoera- 
ture  soundings  and  climatological  salinity  data  (see  Table 
VIII) .  The  Brunt-Vaisala  frequency  is  the  free  oscillation 
frequency   for   a    given   density    gradient    : 

Density  gradients  were  calculated  by  the  equation  of  state 
from  the  observed  temperature  profiles  sampled  by  XBT  data, 
and  assuming  a  constant  vertical  salinity  profile  of  35.6%o 
in  the  deep  water  and  of  35.4%c  on  the  continental  slope 
from  climatological  data.  Brunt- Vaisala  frequencies  wera 
then  computed  in  the  first  300M  and  eventually  blended  to  a 
deep    mean   climatological   profile    found   in    Gould    [19731. 

The  eigenf unct ions  solving  the  Sturm  Liouville  problem 
of  vertical  modes  give  the  vertical  dependence  of  the 
velocity  fields  for  anv  f T*ecrue.'i r v  of  ^  ","arT,5i~,;,n»ai"tial 
oscillation.      An    infinite   number   of    solutions    exist    yielding 
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vertical  distributions  characterized  by  one  or  more  phase 
reversals  in  horizontal  speeds  with  dap the  (Fig  6.1  and 
6.2)  .  For  eaci  of  these  modes  intarnai  wavelengths  ar.d  phase 
speeds  can   be   found    as  a    function    of    the    forcing    frequency. 

Computations  were  accomplished  using  a  numerical  solu- 
tion [Mooers  1971]  to  the  eigenvalue  problem.  Starting  with 
the  boundary  condition  a4,  the  surface,  successive  correc- 
tions cf  the  eigenvalue  are  made  at  each  iteration  to 
satisfy   the   bottom    boundary    conditio-).. 

The  solutions  are  defined.  »ithin  a  multiplicative 
factor.  An  estimate  of  the  vertical  distribution  of  ampli- 
tudes values  is  then  necessary  to  resolve  this  incertitude. 
A  simpla  best  fit  of  tha  ampiituiB  values  from  available 
data  was  made  to  the  first  baroclinic  mode.  Vertical  ampli- 
tudes at  the  semidiurnal  frequency  were  axtracted  from  the 
current  meter  data  (Cf  IV-C)  :  tamoerature  fluctuations  as 
related  to  tha  mean  vertical  thermal  gradient  at  the  corre- 
sponding depth  gave  an  estimate  of  the  vertical  excursions 
of  the  isotherms  (as  shown  in  IV-Ei  .  (The  variations  of  the 
current  meter  pressure  gauge  showad  that  vertical  oscilla- 
tions  of   thc   array    itself   were   negligible) 

The  second  mode  vertical  profila  is  coherent  with  those 
observations  and  is  consequently  assumed  to  ba  the  mean 
component  of  internal  anargy  at  the  semi-diurnal  frequency 
(Fig.  6.3).  Maximum  vertical  amplitudes  are  inferred  to  be 
a*  2  0  001  and  of  the  order  of  5  0M,  a  level  of  no  vertical 
movement     (medal    node)    existing   at    1  3  3  3  M  . 

The    analysis      of    the    3T      soundings    exhibits      an    offshore 


spacing      of   perturbations 


4 5 KM      corresponding   to      phase 


speeds  of  1 M/3  for  a  semidiurnal  generation  and  wavelengths 
of  7.2KM  in  the  open  sea  (90MIN  period)  from  tha  first  mode 
dispersion    analysis. 


Q? 
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In  conclusion,  the  internal 
mod?,  does  not  appear  d n  the  SAR  image,  neither  in  the 
spatial  analysis  of  the  the rmost rupture.  Nevertheless,  peri- 
odic   disturbances   could    be    generated    at    a    semi-diurnal    rate, 

consisting  of  high  frequency  waves  which  are  discussed  now. 
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C.        HIGH    FREQUENCY    INTERNAL    WAVES 

The  wavelengths  in  ths  internal  wave  packets  are  of  the 
order  of  3KM,  yielding  frequencies  of  the  order  of  70MIN 
which  are  consistent  with  the  internal  wave  climatology  (Cf 
III-F).  The  spacing  of  the  high  frequency  wiv=  packers  as 
seen  in  the  propagation  area  of  the  3AR  image  is  approxi- 
mately   23KM. 

On  the  continental  shelf,  the  periodic  spacing  in  the 
patterns  is  about  28KM,  which  is  close  to  ths  semidiurnal 
internal    wavelength    computed   for    shallow    water. 

1 .      Estimation    of  the    Mixed   Lay_er    Seoth 

It  is  assumed  that  the  clisa tological  iensity  struc- 
ture of  the  Bay  of  Biscay  is  known,  and  that  frequent  XBT 
surveys  as  well  as  remote  sensing  oi  the  sea  surface  temper- 
ature, allow  estimation  at  any  time  of  ths  mean  density  jump 
at  the  thermocline.  Starting  witi  only  this  hypothetical 
element  and  Baking  abstraction  of  the  actual  in-sifu  lata 
available,  it  is  attempted  to  compute  the  thermocline  depth 
as  a  function  of  the  geometric  characteristics  of  internal 
waves    as    derived    from   the    SAR    image. 

Assuming  a  dim  at  3  logical  density  jump  of  1.5x10  3 
at  the  thermocline,  and  observing  separation  between  wave 
packets  of  23KM  in  the  op=n  sea  area  of  the  SAR  image  yields 
an  estimate  of  the  mixed  layer  depth  of  20M.  A  variation  of 
10^  in  the  assumed  density  jump  would  result  in  10*  differ- 
ences in  this  estimate,  whereas  changes  of  13^  in  the  inter- 
packet  distances  would  result  in  2d%  changes  in  this  derived 
mixed  layer  depth.  Consequently,  if  the  absolute  estimate 
of  the  thermocline  depth  is  uncertain,  comparison  between 
interpacket  distances  would  be  a  very  sensitive  way  of  esti- 
mating  time    variations. 


9U 


2»      Estimation    of  the   Internal    jlill§   Amglitqda 

Amplitudes  of  the  nonlinear  short  internal  waves 
were  computed  by  a  method  described  and  applied  on  a  similar 
SAR  image  by  Apel  and  Gonzalez  [1993].  The  firs-  assumption 
made  was  that  in  the  "generation  area"  {continental  slope) , 
the  internal  waves  visible  on  the  SAR  image  could  reasonably 
be  considered  as  solitons.  By  definition,  a  scliton  is  a 
perturbation  conserving  its  shape  and  general  characteris- 
tics (except  phase)  ,  even  after  collision  wi'-h  another  one. 
This  characteristic  is  apparent  in  the  vicinity  of  point  P1 
for  example  (Fig.  5.3)  ,  where  ssveral  wave  packets  are 
superimposed  without    any   notable   mataai    interference. 

This  conservative  nature  rf  solitons  is  linked  to 
the  fact  that  the  linear  phase  spe  =  i  dispersion  is  counter- 
acted by  a  nonlinear  acceleration  dependent  on  the  amplitude 
itself.  This  balance  is  -xpressed  in  the  last  two  terms  of 
the    Kortweg    -   de   7ries  equation: 


"t|  +  co\^r^+^ur-o 


y 


whsre  "^  "  is  the  amplitude  of  the  solitary  wave  propagating 
with  time  "t"  in  the  "x"  direction,  "X"  is  the  coefficient 
representinQ  the  linear  dispersion  and  "o^"  is  the  nonlinear 
coefficient    and    is    given    by 

where  "W"  is  tha  vertical  eigenf unction  or  structure  func- 
tion, "h  ,  "  and  "h^"  are  the  dsp'-hs  of  the  layers  in  a 
simplified    discrete    model. 

From  the  mean  density  stratification  observed  a-1:  the 
time  of  the  SAR  image,  and  using  the  vertical  modal  analysis 
described  above,      linear    piase    spoils    "C   "    for    aach    observed 
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wavelength,  and  a  phase  spaed  limit  "c  ,f  for  infinitely  leng 
internal  waves  in  a  nou  rotating  ocean  (steady-state  solu- 
tion   to   the   Stur  m-L  iouvilia    problem!     were    computed. 

In  ona  givan  packa1:  ,  the  wavelength  stretch,  ("i")r 
was  defined  as  the  ratio  of  the  observed  wavelength  to  -he 
linear  Theoretical  one.  From  the  tha:>ry  [Apel  and  Gonzales 
1933],  this  ratio  is  equal  to  "2K  (a2  ) '* ,  twice  the  complata 
elliptic  integral  of  the  firs-  fcini  with  argument  m2  ;  the 
argument  m2  baina  itself  a  function  of  the  soli-en's  ampli- 
tude   A: 

?or  any  wavelength  it  is  then  possible  to  extract 
tha  noniinearity  dua  to  a  given  amplitude.  Or,  conversely, 
observing  a  nonlinear  dispersion  in  tha  patterns  of  the 
internal  waves  (by  measuring  the  variation  of  wavelength  in 
a  given  pacJcat)  ,  it  is  pcssibii  to  derive  ar.  =s-imaoed 
amplitude   of  tha    waves. 

ti^.i  was  *aken  -.0  ba  of  ths  order  of  5::1Q-2  m  l,  ("«C"  '*»20H>  , 
"C0"    was    found    to   ba    0.55^/S. 

For  each  obsarvad  wavelength  in  a  givan  packet ,  the  corre- 
sponding phase  speed  can  be  computed,  and  the  wavelength 
stretch  can  ba  msasured.  A  value  for  "m2"  is  derived  from 
which   wave   amplitude    is  calculated   is 


k  ^     cq-cp  .     3  m1 


^o  =< 


For  an  average  stretching,  taxar.  froa  the  wavegroups 
visibia  in  "ha  generation  arca,  iaoiitudas  of  15^  for  tha 
laading   waves   were    found. 
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D.  INTERNAL   WA7E    INTERACTION    WITH    5H2LFBR2AK    STUCTURE 

Two  different  kinds  cf  geographic  chiig^s  in  the 
internal  wave  group  characteristics  are  apparent  (Cf  7-3-3). 
The  general  shortening  in  wavelength  will  be  assumed  here  to 
be  either  a  Doppler  shift  due  to  a  current  shear  or  a  phase 
speed  compression  due  to  a  horizontal  density  gradient  (e.g. 
as  in  a  front  or  a  mixing  area) .  Refraction  of  the  crests 
resulting  in  a  deformation  from  their  original  cylindrical 
spreading  as  in  groups  northeast  of  ?1  and  ?2 ,  can,  on  tha 
ether  hand,  be  more  surely  associated  with  a  strong  near 
surface   current.  Estimates   of      -he    direction      of    a      shear 

surface  current  as  given  by  characteristic  refraction 
patterns  are  indicated  on  figure  5.4  rhey  give  some  indica- 
tion of  a  general  southwest  ward  alongshore  flow  to  the  south 
of  the  generation  area  (on  the  continental  slooe)  ,  opposite 
to  the  more  localized  currents  at  che  shelf  break,  and  to 
those   south   of    the    slope    (Fig.    5 .  4»  . 

E.  SWELL    REFRACTION    AND    53ELF3REAK    CURRENT 

A  subjective  analysis  of  tas  swell  refraction  was 
attempted.  No  Fourier  Trar.s  forms  were  made  of  ohe  SAR  imags, 
but  in  nost  parts  of  it,  the  clarity  of  the  swell  field  is 
such  that  visnal  estimation  of  its  direction  and  (with  less 
precision)  cf  its  length  can  be  male.  It  was  thought  inter- 
esting to  try  to  apply  recent  Theories  on  wave/current 
interaction  no  this  particular  case,  where  some  severe 
refraction  is  obvious  and  the  ores?ice  of  an  interacting  jet 
probable.  The      theory      initial!/      developed      by      Phillips 

[198  1],  has  been  applied  to  a  SA5  image  of  the  Gulf  Stream 
by  Meadows  [1933];  it  was  transposed  with  drastic  simplifi- 
cations   to   this    case. 
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In  the  same  way  that  the  Snail's  law  (loi  da  Descartes), 
de  scribes  the  r  a  fraction  3f  a  ray  2  a  countering  a  gradient  of 
propagation  in  a  medium,  the  following  equation  describes 
the  refraction  of  swell  of  a  given  period,  propagating  with 
a  phase  velocity  "c"  through  an  interface  between  two  sea 
surface  areas  with  surfacs  currents  af  "  ua  "  and  ,MA(";  "c<0" 
and  "o(,"  are  the  angle  of  incidenca  aid  of  refraction  of  a 
swall    "ray"    (Fig.      5.5). 

us*.,  «•*  *■• 


Assumina  that  the  cause  of  the  observed  refraction  is  a 
shear  current  directed  along  the  shalfbreak  (or  at  Isast 
parallel  to  the  isobaths)  ,  tha  angla  of  incidence  and  the 
angle  of  refraction  (jus-  after  ths  first  visible  break  in 
the  wave  crests)  were  measured.  ras  currant  spaed  (equiva- 
lently  the  let  velocity)  was  estiaatad  to  be  O.SM/S.  The 
accuracy  obtained  cannot  be  more  than  50%,  but  -.he  basic: 
observation  seems  reliable  enough  ,  at  least,  to  suggest  the 
presence  of  a  jet  and  to  localize  the  position  of  the 
maximum    speed   at   the    continental    braalc. 
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711.    TENT ATI  73    BX^HAPOLAriON    £0    5SHZHAL    OCEANIC    FUTURES 

A.       ITJTSBNAL   5JA7E   GENERATION 

"In   the     coastal",  region      internal    tides      generated   or.     a 

continental  slope  are  standing  waves"  [Rattray  1959].  The 
modal  analysis  relative  to  the  density  structure  observed  at 
tha  time  of  the  SAR  image  shows  wavelengths  of  the  order  of 
27K3  for  the  semidiurnal  first  mode  internal  tide  on  the 
shelf.  This  wavelength  is  comparable  with  the  spacing  of 
patterns  visible  on  the  3AR  image  iD::h  of  -he  shelfbreak. 
It  can  therefore  be  assumed  that  thssa  patterns  correspond 
to  the  interaction  of  another  internal  perturbation  periodi- 
cally modulated  in  phase  with  ridgss  and  troughs  in  an 
internal  tide.  Obviously,  the  assumption  of  a  perfect 
standing  internal  tide  would  only  be  true  if  fche  boundaries 
of   the   shelf  sea   acted  as    perfect    reflec-ors. 

Wavelengths  of  the  first  internal  bar  oolitic  modes  in 
deep  water  ace  170km  for  mode  one,  76KS  for  mode  two  and 
&4K3  for  mode  three.  So  direct  manifestation  at  these 
length  scales  is  visible  on  the  SAR  image.  Nevertheless,  the 
high  frequency  internal  wave  packets  give  an  indirect  demon- 
stration of  the  high  energy  loss  enoountered  by  the  internal 
tide  at  the  irregularities  in  the  shelfbreak.  Located  at  the 
heads  of  the  narrow  canyons,  oscillations  in  the  water 
column  are  estimated  (for  this  specific  case  of  the  spring 
tide)  to  be  of  the  order  of  15M  which  are  as  suited  to  be  lee 
waves  forced  by  the  ebbing  tidal  current.  The  internal  wave 
packets  generated  by  this  mechanism  are  localized  at  the 
canyons  extending  toward  the  top  of  sIod*  to  at  least  the 
10  0  0  M  isobath  and  occur  in  bursts.  They  eventually  develop 
anisctropicaily    due    to   the    presence    of    a    shaifbreak    fro?.-    to 
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the  north.  The  general  direction  of  propagation  is  toward 
the  south-west  perpendicul ar  to  the  continental  slope-  \ 
thorough  discussion  of  t.hs  possible  genera-ion  mechanisms  is 
given  by  Pir.gree  [1981].  The  specific  influence  of  narrow 
canyons  has  been  recently  enligh-enei  by  a  laboratory  exper- 
iment  [Baines    1983].. 

B.  TIDAL  MIXING  ON  THE  CONTINENTAL  SLOPE 

Tidal  mixing  on  the  continental  slope  is  demonstrated  by 
the  divergence  in  isotherms  in  the  juasi-synoptic  X3T  anal- 
ysis, (Fig.  4.10  and  4.13)  .  end  is  in  good  agreement  with 
other  measurements  made  i.i  summer  and  fall  at  the  same  loca- 
tion. In  Pingree  (1981) f  the  various  hypothesis  related  to 
this    mixing   are    developed. 

C.  SHELFBREAK    FRONTS    AND    JETS 

A  new  insight  on  the  sheifbrea'c  circulation  is  given  by 
the  SAR  image.  Parallel  bands  t  f  divergence-convergence 
which  compress  internal  waves,  can  be  associated  with  the 
thermal  structure  inferred  from  3T  lata  and  illustrated  bv 
the  following  simplified  scheme  of  circulation  (Fig.  "7.  1)  . 
Over  the  shelf break ,  the  dynamic  heights  tend  to  be  lower 
than  on  either  side  as  a  consequence  of  the  tidal  mixing. 
Therefore,  a  strong  gecscrophic  flow  extends  to  the  thermo- 
cline,  follows  the  isobaths,  and  is  confined  to  the  shallow 
part    of      the   continental   slope.  Trie    flow      is    consequently 

toward  the  northwest  on  the  uooer  oart  of  the  sicoe  and 
toward  the  southeast  in  the  off  sis  re  half.  The  idealized 
transect  (Fig  7.1)  ,  accounts  for  the  localized  and  strong 
northwest  flow  and  explains,  by  stretching  of  the  thernto- 
ciine  offshore,  the  lower  values  ?stimated  ftr  the  count  er- 
current.  It  is  nevertheless  obvious  that  any  '<ind  of 
analytical      model  oased      :n      geostrophy    (steady-state)         and 
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alongshore  homogeneity  cannot  be  fully  appropriate  w  he- 
applied  to  ths  Celtic  sea  shelfbreate.  In  -his  ar=a  currents 
are      dominant  and     limited    geographically     by      the    coast     of 

3ritanr.y,  bathymatr ic  irregularities  are  essential  in  the 
internal  generation  processes,  and  aoniinear  effects  between 
the  steady-state  circulation  and  the  trial  movements  may  be 
expected. 

The  most  informative  indication,  of  the  presence  cf  a 
strong  alongshore  jet  are  the  swell  refraction  ani  the 
sudden  increase  in  wave  height  (17-8,  7-C,  71-5).  The  coin- 
cidence of  these  phenomena  give  reliable  evidence  of  a  jet 
which  may  be  related  to  the  poleward  eastern  boundary 
current    discussed   by    Swallow   [1977], 
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VIII.     APPLICATION    TO    2NDER£AT2R    ACOOSTICS 

Of  ths  host  of  variables  af farting  sound  propagation  in 
the  ocean,  the  temperature  profile  has  the  strongest  influ- 
ence. Not  only  the  vertical  gradie.it  profiles,  but  also  the 
horizontal  gradient  fields,  determine  the  widely  varying 
shapes  of   the  insoni fication    volumes. 

Hence  -enDorai  and  spatial  t  emperatura  variability 
should  be  taken  into  consideration  wher.  forecasting  any 
acoustic  quantity,  from  one  single  climatic  or  observed 
temperature    profile.  Two    approaches   are    possible      in   -his 

matter:  (1)a  deterministic  approach  would  be  possible  only 
in  the  cases  where  the  water-mass's  vertical  novements  are 
known  as  functions  of  time  and  geogoa paical  locations,  well 
er.ouah  to  allow  oh-  forecast  of  effective  variations  of  the 
acoustic  transmission  loss  with  time,  location  ar.  d  orienta- 
tion of  the  sound  •  beam,  and  (2)  a  probabilistic  approach 
would  attempt  to  model  the  oscillati ons  without  any  phase- 
lock  and  would  only  yield  sstimates  of  the  acoustic  propaga- 
tion   variability. 

Obviously,  in  the  oommon  case  where  the  geophysical 
variables  vary  on  a  short  time  scale  and  without  any  period- 
icity, -.he  probabilistic  approach  is  the  only  one  possible. 
On  the  other  hand,  with  laraer  time  and  soace  scales,  an 
effort  in  sampling  the  actual  geophysical  environment  is 
eminently  beneficial.  In  the  3ulf  3orea.ii  region,  for 
example,  knowledge  of  the  positions  of  cold  and  warm  eddies 
has  proven  instrumental  for  acoustio  forecasts.  In  the  case 
of  shorter  scales,  but  where  oer iodic  motions  are  related  to 
a  known  generation  mechanism  (e.g.,  -ides)  ,  the  same  kind  of 
hind  casting    may    be    possible. 


108 


It  is  attempted  here  to  shew  how  the  ocean  variability 
observed  in  the  Say  of  3iscay  may  affect  sonar  detection 
ranges,  and  to  assess  a  possibility  of  analysing  the  depen- 
dence of  these  ranges  or."  relative  positions  of  scarce  and 
receiver  vs.  the  internal  propagating  wave  field,  as  well  as 
on  the  direction  of-  propagation  relative  to  the  internal 
wave  orientation.  The  internal  variability  taken  as  a  basic- 
assumption  is  an  idealized  extrapolation  of  the  results 
gathered  from  the  study  of  the  SA3  image  on  August  20  1978. 
Therefore,  it  applies  only  to  this  given  particular  situ- 
ation,   and   with    a  wide  range  of  uncertainty. 

This  approach  may  prove  feasible  in  similar  cases  where 
internal  waves  can  be  occasionally  sampled,  either  from 
space  using  a  SA3,  or  by  adequate  in  situ  sampling  techni- 
ques based  on  a  thorough  knowledge  of  the  internal  wave 
climatology. 

The  influence  of  an  internal  wave  field  on  acoustic 
propagation  is  especially  critical  at  high  acoustic  frequen- 
cies, where  diffraction  and  leakage  are  minimal  and  where 
pure  gecmetrio  ray  tracing  can  be  applied.  The  insonified 
volume  will  then  drastically  vary  Ln  the  three  dimensions. 
Long  internal  waves  will  act  as  long  scale  time  variations 
in  an  observed  or  predicted  SS?  (sound  speed  profile).  The 
effect  of  short  waves  can  be  approached  stochastically  and 
lead  to  establishing  statistical  conclusions  on  the  range 
variability.  The  variation  in  propagation  will  be  dramati- 
cally dependent  on  (1)the  relative  position  of  source  (or 
receiver)  with  respect  to  the  ware  field,  (2)  the  relative 
direction  of  sound  oroDaaaticn  and  wave  field  orientation, 
and     (3) the    mean    S3?    and    its    variability    in    space. 

Two  cases  are  to  be  considered  .  (  1)  A  long  (e.g.  tidal) 
internal  wave  field  for  low  frequencies  (long  scales)  to 
evaluate  variation  on  long  ranges  (of  order  20K3)  of  sound 
propagation:      in   this  case    the    Parioolio    Equation     (PI)     .'lodel 
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has  been  used.  (2)  A  short  (as  sss:  oa  the  3  AR  iraace  of  the 
co 1 1 isen tal  sioce)  free  internal  wave  field  a~  high  frequen- 
cies (short  scales)  to  evaluate  in tsraction  with  a  high 
frequency  acoustic  wave.  In  this  oass  the  HZDJSA  model  was 
run. 

A.       IDEALIZED    WAVE    FIELDS 

The  acoustic  model  *as  evaluated  assuming  that  in  a 
vertical  plans  the  sound  speed  profile  varied  from  a  near. 
profile  in  a  periodic  way  (Fie.  3.1)  .  This  assumed  profils 
was  drawn  from  -he  XBT  data,  and  was  -her.  idealized  and 
smoothed  to  present  a  perfect  surface  mixed  layer t  a  sharp 
thermociine      and      a      climat olegical    deep      profile.  For     a 

constant  SS?  two  oonvergenoe  zones  ire  apparent:  one  due  to 
ths  Mediterranean  outflow  at  1100M,  and  another  due  to  the 
pressure   increase      with   depth.  r.hs    sharp      thermociine   was 

made  to  varv  Le  decoh  without  any  deformation  in  its  thick- 
ness or  its  slope  (Fia.  3.1)  ;  the  siape  of  this  depth  varia- 
tion was  mcdelsd  as  a  periodic  car/?  with  wavelengths  of  10 
and    2  0KM. 

Increasing  -his  internal  wavelength  by  a  factor  " kM  is 
equivalent  to  modeling  the  ray  tracing  in  a  plane  oriented 
with    respec-   to    the    wavecrests,    ai    an    angle    "A"    such   that 

sia  (A)  =1/k  . 

It  was  then  possible  to  calculate  transects  of  the  three- 
dimensional  insonif ication  (or  transmission  loss)  in  various 
re  i  a  t  i  ve  d  ir  2  o  t  ions  fro  m  *  he  son  cos  and  d~*aw  hcriz  on  -13  ] 
diagrams    a-1*    one    given  depth     (Fig.    3.3    and    9.4). 
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3.       TRANSMISSION    LOSS   DERI7ED    ?ROM    THS    PE    MODEL 


The  parabolic  equation  (PE)  nodal  developed  by  NOEDA 
[Brock  1978],  modified  by  NUSC  [Lee  and  Bctseas  19321,  and 
adapted  (with  acre  tnodif icationsj  for  use  at  the  Naval 
Postgraduate  School,  gives  a  normal  mode  solution  to  the 
wave  problem  in  the  form  of  an  elliptic  equation.  The  calcu- 
lations made  with  this  model  must  be  considered  with  seme 
caution  as  ^he  initialization  steo  will  affect  the  rear-* 
field  results.  Principal  significance  should,  therefore,  be 
attached  to  PS  model  calculations  it  the  convergence  zones 
(CZ)  ranges.  In  the  case  where  tae  scarce  was  p laced  ever 
an  internal  wave  crest,  no  significant  directional  deoen- 
dence  was  observed  and  the  transmission  loss  was  equivalent 
to  that  found  in  a  constant  SS?  environment.  Conversely,  in 
the  case  where  the  source  was  placed  above  an  internal  wave 
trough,  the  convergence  cones  were  found  to  be  enhanced  by  a 
varying  SSP.  Horizontal  diagrams  of  transmission  loss  were 
drawn  for  two  levels  for  the  receiver  (15  and  45:1)  by  inter- 
polation between  five  vertical  transmission  loss  curves  at 
various  directions  relative  to  the  internal  wave  field  (Fig. 
8. '4)  . 


C.       RAT    TRACING    WITH    SEDUSA 

The  MEDUSA  Range  Dependent  Ray  Model  was  developed  by 
the  Applied  Research  Laboratories  of  the  University  cf  Texas 
at  A  use in  under  NOHDA  sponsorship  lacing  ?Y30.  MEDUSA  uses 
smooth  functions  (splines)  to  interpolate  between  the 
discrete  3S?'s  (given  as  input)  and  solves  the  various  ray 
equations  numerically  for  a  limited  number  of  characteristic 
rays  (eigenrays) .  Eigenrays  are  located  by  interpolation  in 
range  (not  in  depth  as  in  other  ray  tracing  aodels)  between 
rays  bracketing  the  receiver  in  range.  The  scarce  is 
assumed   to    lie    on   the   vertical   axis    tf    a    cviindrical    coordi- 
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sate      svs-:?!,      around     which  the      geophysical   variables      are 
symmetrical.      Eigenrays  ara    located    by   interpolating   between 

rays  bracketing  'the  receiver  in  range  (net  in  depth)  . 
(Hence  any  iirect  modeling  of  horizontal  refraction  is 
excluded).  The  PROFTRP  program  permits  diagnosis  cf  profile 
artifacts  produced  in  the  geophysical  interpolation  routine. 
RAYPLT   produces    ray    path    plots. 

This   ray     tracing  program    was    issd      in      the    same      way   as 
the    ?E      model   to   comoute      acoustic   transmission      in    vertical 


planes   at   various   anales    with 


arnai    wave    ; i ? ^ s   orien- 


tation. Hu4:  here,  the  features  which  were  investigated  wer- 
at  shorter  raiges;  the  direction  dependence  exhibited  in  the 
horizontal  diagrams  (Fig.  8.3)  relates  to  an  eventual 
refraction  under  the  mixed  layer,  thus  cancelling  (for  rays 
travelling  at  an  angle  with  the  internal  wave  crests)  the 
"shadow  cone".  Results  are  only  puraly  qualitative  and  would 
need  to  be  investigated  from  a  transmission  loss  point  of 
view . 

Medusa  was  used  also  t o  input  the  horizontal  sound  SDeed 
orofile  as  observed  in  -he  same  idealized  internal  wave 
field  along  an  horizontal  in  the  direction  of  wave  propaga- 
tion and  at  the  mean  depth  cf  the  thermocline  (Fig.  3.2). 
The  result  of  ray  tracing  in  the  horizontal  plane  will  be 
valid  only  given  the  assumption  taat  the  rays  do  not  make 
any  depth  excursion  from  -.he  choosea  depth  of  the  horizontal 
SSP.      This    aseumotion   can    be      made    only    because    the    vertical 


ray      tracing      shows    that      there 
channel,      so   that    a    significant 


oa: 


a    surface      sound      soeed 


of   the    acoustic   enercrv 


- 1 1 


wii 


orooaa; 


within  the  range   3f  variation 


the  mixed 


layer.    Above  the  the  thermocline   the  sound  speed  will  not 

exhibit   any   horizontal   gradient  and   no 

occur. 


will 
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D.        GENERAL    RESULTS 

Transmission  losses  (TL)  in  an  idealized  internal  wave 
field  were  found  to  be  eguivalant  to  the  TL  relative  to  a 
constant  SS?  environment  if  the  source  was  placed  over  a 
crest.  In  the  case  where  the  source  is  above  2.  -rough,  and 
for  internal  wavelengths  of  10KM,  the  transmission  loss 
decreases  at  the  convergence  zone  by  10dB  for  2."  insonifica- 
tion  perpendicular  to  -he  crests  (Fig.  3.4)  This  anisotropy 
at  long  ranges,  ev en  if  not  representative  of  any  real  ocean 
environment,  is  at  least  an  illustration  "ha"  passive  detec- 
tion in  a  long  internal  wave  field  will  be  Iramaticaliy 
directior.-depen  den t  and  siagest  the  importance  of  ta!<in~ 
this    factor   into   account    in    any   passive    dstastion    scheme. 

Ray  tracing  shows  that  there  is  a  possibility  of  insoni- 
fication  of  the  "shadow  zone"  beneath  -he  theraocline  with  a 
variable  S3?  structure.  This  would  affect  both  passive  an  1 
active  search  (or  escape)  tactical  patterns  a*:  lew  ranges 
(of  the  order  of  4000!*)  .  Ray  tracing  in  the  horizontal 
suggests  the  possibility  of  an  horizontal  duct  (particularly 
in  deep  and  narrow  internal  troughs) ,  but  the  assumptions 
made  here  are  too  drastic  to  allow  any  conclusion  from  the 
elementary  run  made  with  Medusa.  Further  work  in  this  line 
should,  nonetheless,  help  to  tailor  us?  of  any  elementary 
(single  SS?)  ray  trace  model  to  interpret  some  consequences 
of    a    variable  33?   on    transmission    anisotropy. 
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IX.    SAR    AS    AN    INTERNAL    OCEAS    SAMPLING    SYSTEM 

The  SEASAT  experiment  was  the  first  specific  cceano- 
graphic  use  of  a  SAR.from  space.  Assessments  made  a  oriori 
or.  its  monitoring  capabilities  have  been  justified  by  the 
large  number  of  oceanic  internal  processes  it  observed.  In 
less  than  42  hours,  it  covered  an  area  of  100  million  KM2  of 
the  Earth's  surface  (35  of  them  ovsr  the  oceans),  [ Pravdo 
1983 1. 

The  main  specific  advantages  of  =  5  A?  as  an  ocean  remote 
sensor  system  are  its  all-weather  capability  (in  particular 
with  respect  to  cloud  cover,  night  conditions  and  sea- 
state),  its  high  resolution  and  its  sensitivity  to  the 
hydrodynamical    surface  manifestations    of    internal    activity. 

The  cbvicis  difficulty  encountered  with  an  operational 
use  of  such  a  sampling  system  is  the  computation  time 
required  to  process  the  enormous  mount  of  data  (some  106 
pixels  per  second).  In  the  case  of  an  orbiting  station  (as 
proven  by  *he  3IR-A  experiment  on  the  1982  Shuttle 
missions),  the  data  may  be  acquired  and  exploited  on  board 
and/or  kept  in  memory  for  a  further  ground-processing.  But, 
for  an  unmanned  satellite,  data  recording  and  exploitation 
require  that  an  earth  staoion  be  ao  radio  receiving  distance 
from    it. 

A.       GENERAL    REQUIREMENTS    PO  R    A    SAR    OCEAN    MONITORING    SYSTEM 

To      possess    this      particular      capability      to    samplB      the 
internal    ccean,     a   3AR  system   must    have   some    specific 
characteristics. 
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It  must  jiv?  a  continuous  radar  coverage  over  length 
scales  greater  than  the  mesoscale  features  -hat  must  be 
detected;  This  a  3AR  designed  solely  for  sampling  the  swell 
■ay  well  work  07  sampling  10KM  squares  of  ocean  each  1G0KH 
but    will    be   uiabla    to  observe    internal   wave   fields. 

another  quit?  obvious  characteristic  is  related  to  the 
process  cf  creating  a  synthetic  aperture.  The  orbit  charac- 
teristic of  "ha  spaceborta  vector  must  provide  a  relative 
velocity  with  respect  10  the  earth's  surface,  thus  a 
geosynchronous      satellite      cannot      be      used.  Nevertheless 

repeating  orbits  or  quasi-geosynchronous  orbits  (i.e.f 
near-polar  orbits  with  periods  equal  to  the  earth's  rota- 
tion) could  combine  the  convenience  of  a  soecific  coverage 
over  a  given  meridian  area  and  the  necessity  of  a  relative 
soeed  ccmDonent  to  permit  use  of  the  Doppler  integration 
scheme  [Tomiyasu  1982].  ?or  the  investigation,  of  internal 
features,  a  sun-synchron.o  us  orbit  would  no4-  be  optimal 
because  a  given  area  would  always  be  lonitored  at  the  same 
time    each   day. 

The  radar  frequencies  at  which  the  surface  wave  field 
responsible  for  the  Bragg  effect  will  effectively  interact 
with  internal  features  say  vary  from  the  3-bar.d  (a  few 
centimeters  wavelengths)  to  the  L-band  (a  few  decimeters). 
High  frequencies  should  yield  better  wind-speed  sensitivity, 
but  would  increase  the  difficulty  of  inter oretaticn  since 
the  dynamics  of  gravity-capillary  waves  are  more  complicated 
than    those    of   pure    gravity    waves. 

The  viewing  angle  may  itself  vary  from  at  least  15°  to 
600  (ST2-3,  Spaceborne  Imaging  Radar  of  ^h2  coming  Soace 
Shuttle  flights)  ,  will  have  a  capability  to  vary  mechani- 
cally the  viewing  angle  within  tiose  limits).  It  seems 
plausible  that  if  low  grazing  angles  should  give  better 
results  for  imaging  cf  tilted  surfaoe  swell,  high  viewina 
angles  may  be  aore  purely  sensitive  to  hydrodynamical 
straining   by  internal   phenomenon. 
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B.  13    SITU    COMPLEaENTARY    SIMULTANEOUS    PROBING 

To  analyse  characteristics  of  the  internal 
displayed  by  a  5AR,  the  mean  density  structure  must  be 
known.  It  has  been  shown  (Cf  chapter  VI)  that  -.he  mixed 
layer  depth,  the  amplitudes  of  high  frequency  internal 
waves,  and  the  location  of  fronts  oan  be  inferred  from 
consequent  surface  patterns.  It  would  therefore  be  possible 
to  conduct  such  an  analysis  in  real  tine  if  some  X3T 
sampling  was  sade  nearly  synoptical!-/  and  could  be  blended 
with    climatclogical    data. 

Obviously  time  series  of  temperature  structure  would  be 
more  efficient  still  by  giving  synoptic  information  on  the 
frequency  spectrum  and  vertical  nodes  of  the  internal 
motions.  A  renote,  drifting  thermistor  array  as  described  in 
[M03C  1980],  if  launched  in  an  appropriate  position  with 
respect  tc  the  a  priori  internal  wave  field,  could  give  both 
an  absolute  reference  to  radar  raaote  sensing  results  and 
fill  in  the  gaps  between  two  consecutive  passages  of  the 
satellite.  (A  radar  trarspondor  couii  add,  at  lew  costs,  a 
geographic    reference.) 

C.  OPERATIONAL  PROCBSSINS  AND  EXPLOITATION 

The  analysis  of  a  SA3  image  may  b?  lade  on  two  different 
levels  corr  espo  n  di  c  g  to  different  goals,  cine  constraints 
and    processing    facilities    (Fia.    9.1i  . 

A  heuristic,  manual  analysis  , similar  to  the  interpreta- 
tion bv  meteorologists  of  a  visible  or  infra-red  satellite 
image  of  the  atmosphere,  would  probably  be  affective  if 
Other  archived  images  were  available.  X  simple  comparison  of 
similar  features  is  capaole  of  giving  indications  of  the 
relative  properties  cf  a  given  internal  wave  field,  and  of 
the  geographic  displacement  cf  fronts.  Hough  quantitative 
measurements  of      interpacket    separation,         wavelengths,      and 
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Variation  of  wavelength,  could  be  ised  in  simple  anal 7 ti oral 
models  to  compute  features  of  operational  interest  such  as 
internal  wave  amplitudes  and  directions,  and  location  of 
current    jets. 

On  the  other  hand,  for  routine  exploitation  of  SAP. 
images  of  a  given  area,  the  knowledge  of  climatology  could 
probably  help  to  tailor  processing  algorithms,  aventually 
using  the  raw  data  directly.  Once  the  features  of  interest 
are  defined  by  their  approximate  Locations,  dimensions  and 
average      characteristics ,  the     required      data        could     be 

extracted  as  a  subset  and  processed  according  no  specific 
methods.  For  example,  if  a  front  is  known  no  move  parallel 
to  itself  and  to  be  characterized  op.  a  SAP.  image  by  a  band 
of  high  radar  intensity  dae  no  convergent  enhancment  of  the 
surface  waves,  the  mean  radar  intensity  averaged  parallel  to 
the  front  would  be  enough  to  iocane  its  position  (and  inten- 
sity) .  If  the  presence  of  high  fns^uency  internal  waves  is 
assumed,  a  procedure  to  extract  their  mean  direction  and  the 
distribution  of  their  wavelengths  would  be  to  analyse  a  set 
of  pixels  over  a  limited  area  by  a  two-dimensional  Fast 
Fourier   Transform   technique. 
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X.    SUMMARY    AND    RJCOSSSND ATIONS 

The  SAR  image  gathered  by  SEA3AT  over  the  3a7  of  3iscay  on 
August  20  1973  showed  a  number  of  internal  features  which 
were  generally  in  good  agreement  with  analysed  synoptic  and 
climatoicgical  data  and  could  hav=  alternatively  be*n  used 
to  infer  valuable  sea^truth  „in  operational,  real-tiae 
sampling. 

Wave  packets  propagating  offshore  from  the  shelfbreak 
had  an  inter-  packet  sparing  on  taa  order  of  23KM  corre- 
sponding to  a  semi-diurnal  generation  at  the  main  canyons 
indenting  the  shelfbreak  at  low  oarotropic  tide  (maximum 
ebbing  surface  flow)  .  Speeds  computed  from  observed  density 
profiles      were    0.553/S      ♦/-    0.5.  over      the    slope,         which 

compares  well  with  the  wave  speed  of  3.51H/S  calculated  for 
those   packets      from    the      5AR   image.  These    high      frequency 

internal  waves  have  a  period  of  the  order  of  70MIN. 
Assuming  thcss  waves  became  linear  far  enough  from  -heir 
generation  points,  a  mean  depth  of  the  mixed  layer  of  20M 
was  inferred  and  spatial  /ariatic.ns  ascertained  from  clima- 
toiogical   data    and    observed    group   spacing. 

Drastic  spatial  variations  in  tie  wavelength  observed  in 
a  jz.  ven  crrouD  and  between  different  ones  were  ar.alvsed,  and 
assumed  to  be  due  to  current  sheao  and  horizontal  density 
gradients.  They  coincided  with  a  =har?  thermal  front  at  the 
shelfbreak  and  an  alongshore  surface  geostrophic  jet.  A 
line  of  surface  enhancement  of  3rao;j  waves  corresponded  to 
this  front,  which  was  anaiysei  by  quasi-synoptic  X3T 
sampling.  This  front  is  generally  oresent  in  summer  condi- 
tions  when   a   oorresocndinc    upweili.no:    breaks    the    surface. 
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Nonlinear  descriptions  cf  tha  internal  waves  in  each 
internal  wave  packet  yielded  estimates  of  their  amplitudes, 
probably    as    large  as    15H    on    the   continental   slope. 

A  tentative  application  to  underwater  acoustics  shows, 
under  severe  simplifying  assumptions,  that  detection  ranges 
in  such  a  situation  are  essentially  direction-dependent  and 
that  sonar  search  patterns  (or,  conversely,  submarine  escape 
routes)  could  be  greatly  enhanced  if  these  internal  feature- 
could  be  sampled  on  a  regular  basis  and  taken  into  account. 
Adequate  operational  X3T  sampling  could  give  some  in-situ 
information  on  this  variability.  It  is  recommended  -hat,  in 
areas  where  internal  waves  are  known  to  occur,  consecutive 
launches  of  X3T*s  be  made  it  a  time  (or  distance)  interval, 
corresponding  to  half  a  period  (wavelength),  thus  giving  an 
approximation  of  the  internal  wave  amplitude  and  conse- 
quently of  the  variations  in  demotion  ranges  that  can  be 
expected.  Internal  wave  climatology  developed  using  SAH 
images   would  be    helpful   in    designing    a    sampling    plan. 

The  following  ocean  features  of  operational  interest  can 
be    inferred    from    the    SAR    image   studied   here: 

a)  An  area  of  intense  internal  activity  is  clearly  visible 
on  the  continental  slope  in  the  form  of  short  internal  wave 
groups  propagating  from  *the  heads  of  the  ma  joe  canyons  and 
generated  approximately  at  the  time  of  low  tide  at  3rest. 
These  perturbations  probably  correspond  with  a  larger  scale 
widening  (weakening)  of  the  thermocline  and  periodic  mixing 
processes.  .he  possible  operational  consequences  are 
(1)  drastic  changes  in  detection  ranges  due  to  a  varying 
thermal  structure  in  time  and  space  and  a  Drobable  increase 
in  ambient  noise  level  due  to  bottom  stress  and  turbidity 
currents,  and  (2)  sudden  modifications  in  buoyancy  that  may 
affect   underwater   vehicles     (e.g.,    submarines, torpedoes)  . 

b)  A   well-defined    front,       located    it    the    upper    oart    of   the 
shelf  break,        is      apparent    on      the      3  A3      imace   as      a      strong 
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localised  enhancement  of  the  high  frequency  waves.  It  corre- 
sponds to  the  well-known  shelfbreaic  upweiiing  in  the  summer 
in  the  area,  which  has  been  demonstrated  by  synoptic  X3T 
sampling.  Iasonif ication  through  such  a  front  is  essen- 
tially dependent  on  the  relative  positions  of  source  and 
receiver. 

c)  In  the  open  sea,  several  fields  of  internal  waves  are 
organised  in  groups  spreading  from  the  shelf.  The  general 
direction  of  the  crests  is  parallel  to  the  isobaths  and 
creates  anisctropies  in  sound  speed  structure  and  insonifi- 
cation   volume. 

d)  Other  features,  not  considered  here,  have  eventual 
tactical  importance:  (1)the  localisation  of  areas  of  calms 
and  therefore  mirror  sea  surface  for  visual  search  in  kSV 
and  SAB  (Search  and  Rescue)  missions  and  site  selection  for 
weapon  launch,  and  (2)  the  coastal  refraction  of  the  swell 
giving  a  possible  surf  hindcastiag  capability  for  mine- 
sweeping   and  amphibious   landing  operations. 

The  coming  SI3-B  missions  will  cover  the  Say  of  Biscay 
and  consequently  provide  some  information  on  the  time  vari- 
ability of  the  characteristics  appearing  in  this  single  SAP. 
image.  Combined  international  effort  to  increase  synoptic 
in-situ  data  acquisition  at  this  time  should  be  eminently 
ben  3f icial. 

Spaceborne  imaging  radars  will  supply  new  capabilities 
for  understanding  and  perhaps  even  continually  sampling  the 
evolution  of  the  internal  ocean.  These  capabilities  may 
dramatically  change  the  parameters  of  undersea  warfare. 
Svery  effort  should  therefore  be  made  in  allied  navies  to 
understand  and  master  these  consequences,  from  both  scien- 
tific and  operational  points  of  view,  to  make  sure  that  no 
dangerous  and  destabilizing  consequences  could  result  from 
their   universal    exploitation. 
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